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FOREWORD 

This  documant  is  Volume  X of  the  Interim  Report  seriee  for  the  Passive 
Nosetlp  Technology  (PANT)  program.  A summary  of  the  documents  in  this  series 
prepared  to  date  is  as  follows t 

Volume  1 - Program  Overview  (U) 

Volume  XI  - Environment  and  Material  Response  Procedures  for  Nosetlp 
Design  (U) 

Volume  III  - Surface  Roughness  Effects 

Part  I - Experimental  Data 

Part  11  - Roughness  Augmented  Heating  Data  Correlation 

and  Analysis  (U) 

Part  III  - Boundary  Layer  Transition  Data  Correlation  and 
Analysis  (U) 

Volume  XV  - Heat  Transfer  and  Pressure  Distributions  on  Ablated  Shapes 

Part  I - Experimental  Data 

Part  II  - Data  Correlation  and  Analysis 

volume  V •>  Definition  of  Shape  Change  Phenomenology  from  Low  Tempera- 
ture Ablator  Experiments 

Part  I - Experimental  Data,  Series  C (Preliminary  Test 
Series) 

Part  II  - Experimental  Data,  Series  D (Final  Teat  Series) 

Part  III  - Shape  Change  Data  Correlation  and  Analysis 

Volume  VI  - Graphite  Ablation  Data  Correlation  and  Analysis  (U) 

Volume  VII  - Computer  User's  Manual,  Steady-State  Analysis  of  Ablating 
Nosetips  (SAANT)  Program 

Volume  VIII  - Computer  User's  Manual,  Passive  Graphite  Ablating  Nosetlp 
(PAGAN)  Program 

Volume  IX  - Unsteady  Flow  on  Ablated  Nosetlp  Shapes  - PANT  Series  G 
Test  and  Analysis  Report 
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VoluM  X - SuBMary  of  Exporinontal  and  Analytical  Roaulta 

VoluM  XI  - Analyaia  and  Ravlow  of  tha  ABUS  Conbustion  Taat  raollity 
for  High  Prasaura  Hyparthamal  Raantry  Nosatip  Systama 

Tasta 

Voluma  XXI  “ Noaatip  Tranaition  and  Shapa  Changa  Taata  in  tha  AFPDL  50 
MN  HINT  Arc  - Data  Haport 

Voluma  XIII  - An  Bxparinantal  Study  to  Bvaluata  Haat  Tranafar  Rataa  to 
Soallopad  Surfaoaa  - Data  Raport 

Voluma  XXV  - An  BxpariiMntal  Study  to  Bvaluata  tha  Xrragular  Noaatip 
Shapa  Aaglma  - Data  Raport 

Voluma  XV  - Roughnaaa  Xnduoad  Tranaition  Bxparimanta  - Data  Raport 

Thia  raport  waa  praparad  by  Aarotharm  Diviaion/Aourax  Corporation  undar 
Contract  P04701-71-C-0027.  Volumaa  I through  IX  eovarad  PANT  activitlaa  from 
April  1971  through  April  1973.  Volumaa  X through  XV  rapraaant  contract  afforta 
from  May  1973  to  Dacambar  1974.  Voluma  X aummariaaa  tha  raapaotiva  taat  pro- 
grama  and  daacribaa  improvamanta  in  noaatip  analyaia  oapabllitlaa.  Voluma  XI 
praaanta  an  evaluation  of  tha  ABRB8  taat  facility  in  tarma  of  performing  ther- 
moatructural  and  reentry  flight  aimulation  taating.  Volumaa  XII  through  XV  are 
data  raporta  which  aummariia  tha  axparimanta  performed  for  tha  purpoaa  of  de- 
fining tha  irregular  flight  regime.  Tha  analyaia  of  thaaa  data  are  praaantad 
in  Voluma  X. 
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ABSTRACT 

\ 

Results  from  «xp«rlm«ntel  and  analytical  efforta  parformad  on  tha 
Paasiva  Neiatip  Tachnolo^y  (PANT)  program  batwaan  May  1973  and  Daoambar  1974 
ara  raportadi  Pour  taat  programs  wara  conducted  oonaiatlng  of  one  test 
sarios  In  the  AFFDL  50  MN  arc  haatar,  and  thrae  teat  sarias  In  the  Naval 
Ordnance  Laboratory!  Wind  Tunnel  No.  8.  Tost  objectives  and  nominal  results 
ara  summarised  for  each  test  sarias.  Analysis  of  heat  transfer  and  bound- 
ary layer  transition  data  ara  described  and  changes  to  nosatip  analysis 
procedures  ara  raconmandad.  The  affect  of  modeling  changes  on  nosetlp 
response  predictions  ara  illustrated  for  several  flight  environment  cases. 
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SECTION  1 


IMTEODOCTIOM 

During  th«  period  April  1971  through  April  1973,  the  Aorotherm 
Division  of  Aourax  Corporation  waa  involved  in  tha  development  of  reentry- 
vehicle,  paseiva  noeatip  deaign  tools  under  tha  SAMSO  sponsored  Passive 
Nosatip  Technology  (PANT)  program.  The  effort  included  experimental  and 
analytical  investigations  of  the  fundamental  phenomena  controlling  the 
response  of  nosatip  materials.  Results  of  this  effort  were  reported  in 
tha  PANT  program.  Interim  Final  Report,  dated  January  1974,  Reference  1-1. 
The  signifioant  findings  consisted  of  the  following! 

e Tha  importanoe  of  the  nosatip  surface  roughness  on  boundary 
layer  transition  and  heat  transfer  was  demonstrated  through 
wind  tunnel  simulation  experiments. 

e Hind  tunnel  simulation  experiments  also  demonstrated  that  the 
surface  temperature  to  boundary  layer  edge  temperature  ratio 
significantly  influences  the  location  and  ocourrenoe  of  bound- 
ary layer  transition. 

e Nosetip  ablation  tests  in  both  wind  tunnel  and  hyperthermal 
facilities  indicated  that  surface  roughness  patterns  (l.e., 
scallops)  which  affect  heat  transfer,  develop  during  turbu- 
lent ablation. 

a Consideration  of  roughness  effects  on  heat  transfer  indicated 
that  the  micromechanical  particulate  component  of  graphitic 
nosetip  mass  loss  is  negligible  and  that  graphite  mass  loss 
and  surface  temperature  data  can  be  predicted  utilizing  equi- 
librium thermochemistry  theory. 

e Nosetip  shape-change  data,  primarily  from  wind  tunnel  tests, 
indicated  that  nosetip  ablation  can  produce  irregular  (slen- 
der or  asymmetric)  shapes  which  are  prone  to  gross  failures. 

Computational  procedures  developed  to  model  the  above  findings  were  only 
partially  successful  due  in  part  to  the  tendency  of  nosetips  to  develop 
an  irregular  shape  in  seme  environmental  regimes.  Further  efforts  to 
better  understand  boundary  layer  transition,  roughness  pattern  effects 
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on  hMtlng>  «nd  irrogular  sh«p«  onsot  w«r«,  tharaforo,  raoonaandad  to 
8AN80.  Thaaa  additional  aotivltlas  wara  par£ormad  undar  PANT  program 
work  statamant,  paragraph  4.2.5,  batwaan  Hay  1973  and  Dacambar  1974  and 
ara  doeumantad  in  Intarim  Raport  Volumaa  XX  through  XV. 

Tha  ganaral  objaotiva  of  Task  4.2.5  was  to  dsflna  tha  raantry 
ragliiios  which  produea  irragular  nosatip  shapas.  In  order  to  satisfy 
this  objaetive  a siaraa  of  wind  tunnal  and  high-prassura  ablation  tests 
was  ragulrad.  Tha  tasts  ware  designed  to  obtain  oalorimater  and  shape- 
change  data  for  geometries  and  envionmants  of  interest.  The  specific 
objectives  of  the  four-test  series  ara  given  below i 

50  m Arc  Heater  Teats  - Test  ATJ-8  graphite,  Mod-lll  carbon/ 
carbon,  and  reverse  chevron  carbon  phenolic  nosatip  models 
in  the  50  MH  arc  using  a * 3 nossle  to  assess  tha  transi- 
tion behavior  and  shape-change  response  of  "real*  materials 
in  a hyperthemal  environmant. 

Series  H Hind  Tunnel  Tests  - Oenerate  heat  transfer  data  for  a 
calorimeter  replica  of  an  LTA  scallop  roughened  surface  for 
comparison  with  inferred  heat  transfer  data  and  predictions. 

Series  I Wind  Tunnel  Tests  Generate  la#  temperature  ablator 
(liTA)  shape  change  data  to  determine  the  effect  of  relevant 
nosatip  parameters  on  the  formation  and  extent  of  irregular 

shapes . 

Series  J Hind  Tunnel  Tests  - Generate  heat  transfer  data  which 
show  the  effect  of  surface  roughness  and  modal  size  on  bound- 
ary layer  transition. 

The  50  MW  tests  were  conducted  in  the  RENT  leg  of  the  50  MH  arc 
facility  at  the  Air  Force  Materials  Plight  Dynamics  Ledioratory.  Test 
Series  H and  J were  conducted  in  Tunnel  No.  8 at  the  Naval  Ordnance 
Laboratory  (NOD  from  February  21  to  February  26,  1974.  Test  Series  I 
was  conducted  at  NOL  from  March  27  to  April  4,  1974. 

These  teat  series  complement  Series  B and  G which  ware  directed 
at  obtaining  pressure  data  on  nosetlps  of  irregular  shape.  A brief 
description  of  these  two  test  series  is  given  below i 

Series  E (Task  4.2.3)  — Generate  high-frequency  pressure,  acceler- 
ation, and  shock-shape  stability  data  to  establish  tha  occur- 
rence of  high-frequency  flow  field  pulsations  on  ablatsd  nose- 
tip  shapes. 
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SariM  0 (Task  4.2.7)  - Osnsrats  high-fraquanoy  praasura,  acoalara- 
tion,  and  ahook-ahapa  stability  data  to  datsmilna  tha  snvalopa 
of  raalialtio  nosatlp  shapaa  which  indues  high-fraquanoy  flow 
fiald  pulsations. 

Tast  Sarias  E and  Q wars  oonduetad  in  NOL  Tunnal  No.  8j  Sarias  E 
was  eonduotad  In  Ootobar  of  1972  and  Sarias  6 was  oonduotad  from  July  25 
to  August  1«  1973.  Sarias  E and  0 data  ara  raportad  in  Eafaranoa  1-2. 

Volums  XI  of  tha  Xntarlm  Raport  sarias  dooumants  an  ancillary  analy- 
sis of  calibration  data  and  raantry  simulation  oapabilitias  of  tha  ABRES 
Combustion  Tast  Facility  at  tha  Air  Foroa  Rockat  Propulsion  Laboratory  and 
is  not  suimnarisad  harain.  Data  raports  for  tha  four  tast  sarias  parformad 
undar  Task  4.2.5  ara  Inoludad  as  Volumas  XXI,  XXIX,  XXV,  and  XV.  Xn  this  vol- 
ums, tha  rasults  of  tha  various  nosatip  axparlnants  and  tha  associatad 
improvamants  in  rasponaa  modaling  procaduras  ara  summarisad.  in  Saotion  2, 
tha  50  MW  arc  ablation  tast  sarias  la  sunmarizad.  In  Saotion  3,  tha 
Sarias  H wind  tunnal  tast  rasults  and  othar  scallop  offsets  data  ara  pra- 
santad.  The  Sarias  X shapa-ehanga  data  and  rasponsa  trands  ara  suirnnarizad 
in  Saotion  4,  and  tha  Sarias  J data  ara  raviawad  in  Saotion  5.  Saotion  6 
ineludao  a rsviaw  of  rasponsa  modaling  improvamants  and  indioatao  tha 
ovarall  affaots  of  tha  improvamants  on  tha  pradlctions  for  savan  basalina 
flight  easas.  Conclusions  and  raoommandations  from  the  program  ara  ra- 
viawad in  Saotion  7. 
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SECTION  2 


AFFDL  50  MM  ARC  TESTS 

The  Fall  1973  PANT  teat  sariaa  at  tha  AFFDL  50  MM  RENT  facility 
was  oonduetad  for  two  purposes t 

1.  Nosatip  shape  change  tasting  in  a ragims  predicted  to  prn- 
duoa  irregular  shapes 

2.  Tasting  to  obtain  data  on  tha  affaets  of  nose  radius  and 
aatarial  on  transition. 

Tha  tast  sariaa  included  exposure  of  16  models  during  six  arc  runs  plus 
exposure  of  one  additional  model  during  a stdbssguant  test  series.  The 
teat  oiMrationa  ware  generally  smooth  and  data  aequisition  goals  wars  met 
although  minor  problems  in  model  design  were  discovered.  Details  of 
the  test  plan  and  data  are  documented  in  References  2-1  and  Volume  xil 
of  this  series  of  PANT  reports,  respectively.  Section  2.1  gives  a brief 
description  of  tbs  testing,  and  Section  2.2  discusses  the  nosatip  rs- 
sponse  results. 

2.1  TEST  SUMMARY 

Irregular  nosatip  shapes  have  been  observed  in  numerous  PANT  wind 
tunnel  tests  of  low  temperature  ablator  models  (see,  for  example.  Refer- 
ences 2-2,  2-3,  or  Volume  XIV).  The  formation  of  these  undesirable  shapes 
appears  to  requira  at  least  the  following  two  conditions t 

e Boundary  layer  transition  should  occur  forward  of  the  sonic 
point  on  the  nosetip 

e The  flow  field  Mach  number  must  be  high  enough  to  produce  a 
strong  curved  shock  which  can  influence  heat  and  mass  transfer 
rates  and  nosetip  shape  change. 

As  a result  of  the  emphasis  in  most  hyperthermal  materials  tests  on  maxi- 
mising impact  pressure,  flow  field  Mach  number  is  generally  near  one,  and 
the  second  requirement  is  rarely  net.  Thus,  development  of  irregular 
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•hap«i  ha*  not  boon  aaan  on  ground  toata  of  graphltio  matarlala.  daaod 
on  analyaoa  prasantod  in  Rafaranoa  2-4,  it  was  datarmlnad  that  taatlng 
in  tha  MKMT  facility  using  tha  0.9-inoh  throat/2. 39- inoh  axit  diamatar 
oontourad  notala  would  bo,  at  most,  marginally  capabla  of  producing  tha 
conditions  likaly  to  form  irragular  ahapas. 

Tha  other  purpose  of  tha  tasting  was  to  provide  additional  data 
on  material  and  noao  radius  affects  on  transition  in  a hypartharmal 
environment.  Tha  PkNT  tranaitKMi  correlation  (Rafaraneo  2-5)  was  de- 
veloped principally  based  on  wind  tunnel  tasting,  and  its  application  to 
hypartharmal  anvironmants  had  bean  axaminad  only  for  graphite  nosatipa 
with  0. 25-inch  nose  radii  (Rafaranoa  2-6).  This  tast  series  was  dssignod 
to  supply  pertinent  data  on  transition  for  oarbon  phenolic,  oarbon/oarbon, 
and  graphite  matarlala  over  a range  of  nose  radii.  Tha  materials  ware 
chosen  to  provide  a variation  of  surface  roughness  charactaristios . For 
transition  tests  of  tha  smaller  models,  higher  stagnation  pressure  con- 
ditions were  necessary.  Thus,  two  of  tha  arc  runs  ware  performed  using 
the  0.9-inoh  throat/1. 49-inch  exit  diamatar  nossla. 

Tha  tast  matrix  as  actually  run  is  shown  in  Table  2-1  although 
data  from  tha  subsequent  tasting  (RTN  52-10)  are  not  included.  Runs  4, 

5,  and  7 ware  directed  towards  tha  shape  change  objective.  Runs  6,  I, 
and  9 primarily  provided  transition  information.  All  models  oxcapt 
three  had  spherical  nose  radii  on  either  a cylinder  or  an  8*  aft  cone. 

Nose  radii  varied  from  0.2  to  0.75  inoh.  Two  models  ware  30*/8*  sharp 
bioonios,  and  one  model  was  a flat-ended,  0.3-lnoh  radius  cylinder. 

In  all  runs  the  bulk  enthalpy  was  lower  than  originally  anticipated 
due  to  arc  heater  modifications  implemented  just  prior  to  testing.  These 
modifications  improved  cooling  and  reduced  erosion  of  the  front  electrode. 
Tha  first  model  (ATJ-8)  tested  failed  oatastrophloally  at  9 seconds  Into 
the  run.  Evaluation  of  the  model  design  indioated  that  epoxy  adhesive 
used  to  mate  tha  ATJ-8  nosetip  to  the  sllloa  phenolic  skirt  probably 
volatilised  and  fractured  the  graphite.  The  testing  sequence  was  re- 
ordered while  vent  holes  were  drilled  in  the  remaining  graphite  models 
of  the  same  design  and  while  Southern  Research  Institute  fabricated  a 
solid  graphite  backup  model  (RTN  52-007,  Sting  3).  All  other  aspects 
of  testing  were  successful. 

2.2  NOSETIP  RESPONSE  RESULTS 

The  results  of  the  shape  change  tests  (Runs  4,  5,  and  7)  provide 
detailed  large  scale  information  on  shape  behavior  in  realistic 
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•nvicoAMmt«  for  oarbon  phonolio  and  graphita.  Tha  modal  bahavlor  wan 
ganarally  Indloatlva  of  fully  turbulant  eondltlona  without  major  aaym- 
aatrlaa  or  irragular  ahapaa.  Thin  waa  partloularly  trua  of  tha  earbon 
phanolio  modala  which,  daapita  attampta  to  produoa  tranaltional  flow, 
alwaya  aahibitad  wall>bahavad  fully  turbulant  ahapaa.  Nlth  tha  ATJ-8 
modala,  amall  protrualona  or  knobs  davalopad  at  tha  atagnatlon  point 
and  in  faot  appaarad  quits  similar  to  pravious  low  tamparatura  ablator 
Bhapa  profilaa  as  illustratad  in  rigura  2>1.  Tha  flgurs  oomparas  an 
ATJ-8  modal  with  a camphor  modal  taatad  at  Ra.  - S.10*/ft  and  - 5 in 
tha  NOL  hyparsonlo  wind  Tunnal  Numbar  8. 


(a)  Run  007,  Sting  2 - Graphita  (b)  PANT  Run  204  - Camphor 

low  tamparatura  ablator 

Figura  2-1.  Comparison  of  Nosatip  Shapa  Profilas 

Both  tha  ATJ-S  and  tha  carbon  phanollc  aphoro/oylindar  modala 
axhlbltsd  two  distinct  racassion  raglmsa.  Aftsr  transition  ocourrad 
at  about  tha  sonic  point,  tha  ahapaa  procaadad  to  sharpan  to  a conic 
configuration  without  substantial  stagnation  point  racassion.  Aftsr 
achiaving  tha  sharpanad  shapa,  racassion  inoraasad  substantially  and 
tha  shapa  ramainad  ralatlvaly  uniform.  Tha  prssharpanad  <30*)  configu- 
rations for  both  matarials  bagan  to  racada  inmadiatsly  without  having  to 
sharpan  thus  achiaving  substantially  graatar  total  racassion  in  a shorter 
time. 

Tha  transition  tests  (Runs  6,  8,  and  9)  produced  satisfactory 
results  axcapt  for  the  two  carbon  phenolic  models.  Tha  carbon  phenolic 
at  tha  low  impact  pressure  appaarad  to  develop  a significant  char  layer 

t and  than  dalaminata.  Tha  carbon  phanolio  models  transitioned  close  to  tha 
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stagnation  point#  but  tha  shapa  ohanga  and  roughnass  wara  doodnatad  by 
tha  dalsmination  oharaotaristioa  rathar  than  tha  fluid  naohanios. 

Tha  raattin^jig  nodala  all  praaantad  usaful  information  on  aisa# 
shapa#  and  aatar\al  propartiaa  affaots  on  transition.  All  modals  tastad 
in  tha  small  (1.^9  inch)  nossla  with  tha  axoaption  of  tha  flat  faoad 
ATJ-8  nodal  axhibitad  immsdiata  trsnaition  which  saamad  to  ramain  for 
a long  tims  in  an  " intar. "widl>\ta  shapa"  pattam  consisting  of  a laminar 
region  axtanding  up  to  ab>?ut:  tha  sonic,  point  and  than  an  obviously  tur- 
bulent region  with  a fine  scala  scallop  pattern.  This  ahapa  a lowly 
sharpens  to  a configuration  which  is  similar  to  low  tamparatura  ablation 
nosetipa#  as  illustrated  in  Figure  2-2. 


(a)  Run  008#  Sting  2 > Graphite  (b)  FANT  Run  216  - Camphor 

low  tamparatura  ablator 

Figure  2-2.  Comparison  of  Transitional  Shapes 

The  large  (0.75  inch  nose  radius)  Mod-IIZ  model  on  Sting  3 of 
Run  006  appeared  to  be  transitional  at  completion  of  testing.  Large 
gouges  ware  evident  starting  at  about  the  45*  location  around  tha  nose. 
Thus#  this  large  nose  radius  Mod-Ill  appeared  to  be  experiencing  incip- 
ient transition  even  at  this  low  pressure  (pt2  “ 12  stm)  teat  condition. 

The  flat  faced  ATJ-S  model  on  Sting  1 of  Run  009  appaarad  to  main- 
tain a laminar  boundary  layer  through  most  of  its  exposure.  This  modal 
experlanoad  vary  little  recession,  in  particular#  during  the  four  sec- 
onds of  exposure#  this  model  did  not  avan  ablate  to  a rounded  blunt 
"laminar"  shape  as  might  be  expected.  It  ie  uncertain  whether  or  not  the 
irregular  ablation  at  tha  cornara  of  this  modal  near  tha  end  of  tha 
exposure  Is  indiostlve  of  transition. 
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As  indlostsd  sbevsf  ths  shsps  ohsngs  rsaponsss  of  tho  nosstlp 
medslt  ttstsd  in  this  50  MM  are  tost  sarlas  are,  in  ganaral,  oonsis- 
tant  with  rasponsa  trands  obaarvad  in  low  tamparatura  ablator  axpari- 
aanta.  It  is  unfortunata  that  tha  fonaation  of  axtramaly  irragular 
(slandar  or  asyonatrio)  nosatip  ahapaa  was  not  obaarvad.  Tha  lowar- 
than-antioipatad  anthalpy  lavala  aohiavad  in  thasa  are  tasts  undoubtadly 
eontrlbutad  to  aissing  tha  irragular  shapa  ragina  in  thaaa  toats. 

Although  tha  transition  data  ara  adaguata  to  yiald  useful  quan- 
titative information  on  transition  of  real  materials,  thaaa  data  have 
yet  to  ba  analysed.  Comparison  between  transition  in  thasa  contoured 
nossla  tasts  with  transition  observed  in  rampad-anvironnant,  flared 
noBSla  tasts  would  provide  an  indication  of  flow  affects  on  boundary 
layer  transition. 
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SECTION  3 

SCAIMP  BPPBCTS  ON  NEAT  TRANSFER 

On*  of  th*  Important  raaults  of  th*  initial  affort*  oonduetad 
under  th*  PANT  program  «a*  th*  discovery  that  surface  scallops  signifi- 
cantly increase  graphite  recession  rat*  above  smooth  wall  predicted 
values  (Reference  3-1) . In  fact  graphite  ablation  data  from  such  divers* 
test  facilities  as  the  APFDL  50  HN  arc  heater  and  th*  ABDC  Rang*  0 
ballistic  range  are  consistently  predicted  only  if  proper  account  of 
scallop  effects  on  heat  transfer  are  included. 

Low  temperature  ablator  (LTA)  wind  tunnel  test  models  conducted 
under  th*  PANT  program  experienced  th*  formation  of  scallops  which  are 
geometrically  similar  to  those  which  form  on  graphite.  An  apparent  in- 
crease in  recession  rat*  was  also  observed  in  th*  LTA  tests.  Based  on 
these  results  a task  was  defined  to  accomplish  th*  following  objectives i 

* Measure  scallop  dimensions  on  LTA  models. 

e Quantitatively  infer  th*  heating  rat*  to  a scalloped  LTA 
surface. 

e Measure  the  heating  rate  to  a calorimeter  replica  of  the 
scalloped  LTA  surface  (Series  H test  series) . 

LTA  models  trere  used  for  these  measurements  as  opposed  to  graphite  models 
for  th*  following  reasons i 

* LTA  models  were  tested  at  hypersonic  conditions  (M.  - 5)  which 
Is  a much  better  flight  simulation  than  most  graphite  ground 
tests  (M^  - 1.7)  (anticipated  to  Influence  cone  scallops). 

* LTA  scallop  sises  are  larger  and,  hence,  can  be  more  accurately 
measured. 

e LTA  thermochemistry  is  simple  (sublimation  - no  chemical  re- 
actions) and  well  defined.  This  aids  in  inferring  heating 
rat*. 

* Test  conditions  are  more  well  defined  (wind  tunnel  versus 
are  heater) . 
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* LTA  nodal*  *r*  largar  than  typical  graphite  model*  (2.S  ver*u* 

0.25  inoh  no*e  radiu*) . Thi*  aimplifie*  the  oonatruotion  of 
the  replica  calorimeter. 

In  order  to  eatiafy  the  atated  objective*  within  the  aeepe  of  the 
program  a eingle  LTA  model  waa  aelected  for  the  majority  of  the  detailed 
meaeurementa  and  teat*.  The  model  waa  aelected  baaed  on  ita  compatibility 
with  the  taak  objeotivea  and  involved  oonaideration  of 

e Hodel  aymtetry 

e Shape  atability 

e Variability  of  raeeaaion  rat* 

* Ability  to  duplicate  LTA  teet  condition  for  a calorimeter 
replica 

Baaed  on  thee*  oonaideration*  BANT  Seriea  D LTA  Run  201  (Reference  3-2) 
waa  aelected.  The  reeult*  of  the  roughnee*  meaaurementa  and  teata  are 
diaouaaed  in  the  following  aeetiona.  Section  3.1  deaeriboa  moaauramenta 
of  acallop  roughnea*  and  inferred  heat  tranafer  rate.  Section  3.2  pra- 
aenta  reaulta  from  the  acallop  replica  calorimeter  teata.  Compariaona 
betvieen  inferred,  acallop  calorimeter,  and  other  rough  wall  calorimeter 
data  are  covered  in  Section  3.3. 

3.1  MEASUREMENTS  FROM  LTA  TESTS 

Two  types  of  measurements  were  made  from  the  LTA  test  datat 
(a)  scallop  roughness,  and  (b)  inferred  heat  transfer  rate.  These  are 
discussed  in  the  following  two  sections. 

3.1.1  Scallop  Roughness  Measurements 

Measurements  of  scallop  roughness  were  taken  from  post  test 
replicas  of  LTA  models  for  the  following  reasons i 

* To  relate  scallop  roughness  heating  augmentation  to  sandgrain 
type  augmentation 

e To  observe  scallop  roughness  trends  with  test  condition 

Figure  3-1  shows  photographs  of  a turbulent  ablating  LTA  surface  from 
the  PANT  Series  D tests  and  illustrates  the  irregular  character  of 
scallops.  Close  examination  of  replicas  of  scalloped  surfaces  shows  that 
even  at  a constant  streamwia*  location  there  is  considerable  variation 
from  on*  scallop  to  the  next,  'ieucuse  of  this  basic  variability  in 

scallop  aise,  conclusions  about  typical  scallop  roughness  must  be  drawn  .j 
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from  • ■tatlstlcally  slgnlficaiit  data  baaa.  Xn  ordar  to  aatabllah  auch 
a data  baaa  for  aaoh  modal  tha  modal  was  dlvidad  into  t«ro  ragionsr  aa 
follows i 

a Tha  high  angla  aubaonie/transonlo  noaa  ragion 

a Tha  low  angla  suparsonio  oona  ragion 

The  two  typos  of  scallop  maasuramonta  ara  dlseuaaad  In  tha  following 
paragraphs . 

Nosa  Region  8oallopa 

Seallopa  In  tha  ausa  ragion  ara  much  amallor  than  eona  aoallopa 
and,  hanca,  greater  oara  must  be  taken  In  thalr  maaauramant.  The  taoh- 
nlqua  used  was  to  aaotlon  plaatar  eaat  raplloas  of  tha  modal.  Those 
raplloas  ware  made  from  RTV  molds  produoad  dlraotly  after  aaoh  LTA  teat 
and  provide  an  aoourata  record  of  surface  details.  Tha  plaster  sections 
wara  than  enlarged  lOX  on  an  optical  comparator,  and  profile  traoas  of 
tha  scalloped  surface  wara  made.  A datum  aurfaea  was  defined  for  aaoh 
profile,  and  roughness  heights  wara  measured.  Figure  3-2  presents  a 
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Figure  3-2.  Typical  Nosa  Ragion  Profile  Trace  (Run  201) 
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typioBl  profile  trace  and  roughneaa  neaeurenente.  In  the  noee  region, 
only  longitudinal  (l.e.,  parallel  to  the  model  centerline)  cute  were 
made.  Thle  la  beoauae  outa  in  the  lateral  direction  (l.e.,  perpendicular 
to  the  centerline)  interaect  the  aurfaoe  at  a aignlfioant  angle  (approxi- 
mately 45*)  and,  hence,  in  order  to  determine  normal  roughneaa  height 
from  the  profile  trace  a correction  muat  be  made.  Unfortunately,  the 
amount  of  the  correction  in  a function  of  )»oth  the  aurfaoe  interaection 
angle  and  the  roughneaa  geometry.  Xiao,  the  correction  beoomea  inoreaa- 
ingly  larger  aa  the  aurfaoe  Internaetidn  angle  in  increaaad. 

Nona  region  acallop  maaaurementn  were  only  made  for  the  baaeline 
caae.  Run  201.  Four  longitudinal  cuta  were  made.  The  roughnaan  maanura- 
mentn  from  all  longitudinal  cute  were  connidered  together  when  calculating 
the  mean  and  atandard  deviation.  Table  3-1  prenentn  theae  reaults. 

TXBLE  3-1 

NOSE  REGION  SCALLOP  ROUGHNESS  MEASUREHBNTS , 

RUN  201 


Number  of 
Measurements 

Mean  Pealc  to 
Valley  Roughness 
(mils) 

Standard 

Deviation 

(mils) 

31 

27.3 

11.5 

These  roughnesaas  will  be  used  later  when  comparing  acallop  and  sand- 
grain  roughness  augmentation  to  heat  transfer. 


Cone  Region  Scallops 


The  cone  region  was  defined  as  that  region  aft  of  where  the  sur- 
face angle  la  less  than  20*.  The  technique  used  to  measure  scallops 
was  basically  the  same  as  that  used  on  the  nose  except  that  the  scallops 
on  the  cone  are  large  enough  that  the  profile  traces  can  be  made  directly 
from  the  plaster  sections i no  magnification  is  necessary.  In  addition, 
lateral  cuts  of  the  cone  surface  can  be  used  for  scallop  measurements 
because  the  intersection  angle  correction  is  small.  Figures  3-3a  and 
3-3b  show  typical  longitudinal  and  lateral  cone  profile  traces  and  scallop 
measurements. 

Cone  scallop  measurements  were  made  on  four  models  covering  a 
range  of  parameters.  Table  3-2  presents  the  range  of  conditions  con- 
sidered. For  each  model  a minimum  of  four  longitudinal  and  three  lateral 
profile  cuts  were  considered.  The  mean  scallop  dimension  measurement  was 
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rlgur*  3-3.  Typical  Cona  Region  Profile  Traces 
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TABU  3-2 

CONE  raSGXON  SCALLOP  NSA8URSNBNT  MATRIX 


Run* 

No. 

Teat  Conditions 

^o 

(psia) 

^o 

CF) 

H 

• 

(loVft) 

201 

480 

585 

] 

5 

10.02 

203 

775 

975 

10.15 

204 

395 

975 

5.23 

207 

786 

989 

■ 

10.16 

Nos* 

Radius 

(in.) 


Daaoription 


2.5 

2.5 

2.5 

1.5 


Baaslina  modal 

rKfact  of  varying  pras- 
sura  at  oonstant  Rey- 
nolds number 

Effect  of  Reynolds 
number 

Effect  of  nose  radius 


*A11  tests  were  performed  on  PANT  test  Series  D. 


calculated  by  lumping  all  longitudinal  and  lateral  cut  data.  These  re- 
sults are  displayed  in  Figure  3-4  whioh  shows  a comparison  of  the  longi- 
tudinal and  lateral  out  measurements i also,  a net  mean  calculation  using 
both  longitudinal  and  lateral  measurements  is  shown.  The  comparison 
between  the  longitudinal  and  lateral  measurements  is  reasonable  consider- 
ing the  basic  variability  in  scallop  sixes.  Figure  3-4  demonstrates  the 
validity  of  the  assumption  that  the  scallops  on  the  cone  form  a single 
statistical  group. 

As  an  indication  of  the  variation  of  cone  scallop  depth  with  test 
condition,  Figure  3-5  shows  the  mean  scallop  dimensions  plotted  versus 
two  test  parameters.  Figure  3-5a  shows  a plot  of  mean  peak  to  valley 
scallop  dimension  versus  tunnel  supply  pressure.  Scallop  analysis  re- 
ported in  Reference  3-3  indicates  that  if  surface  deformation  is  a pre- 
dominant mechanism  In  scallop  formation  that  cone  scallop  depth  should 
vary  as  cone  pressure.  For  comparison,  a plot  of  scallop  depth  versus 
free  stream  unit  Reynolds  number  is  shown  in  Figure  3-5b.  Although 
more  data  are  needed  to  be  conclusive,  these  two  plots  seem  to  indicate 
that  scallop  depths  correlate  better  with  Reynolds  number  than  pressure. 

3.1.2  Inferred  Heating  Data 

In  this  program,  a method  was  developed  for  inferring  the  heating 
rate  to  an  ablating  camphor  surface.  Basically,  the  technique  involves 
combining  measurements  of  surface  recession  rate  with  material 
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rigura  3-5.  Variation  of  Paak  to  Vallay  Bcallop  Dimanaion 
With  Taat  Condition 
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tharmoehwklstry  and  oaloulatad  ad^a  oonditiona  to  infar  tha  haating  rata 
raquirad  to  match  tha  obaarvad  raoaaalon  rata.  It  ahould  ba  amphailsad 
that,  for  camphor,  raoaaalon  rata  la  diractly  proportional  to  haating  and, 
tharafora,  tha  accuracy  of  tha  Infarrad  haating  la  proportional  to  tha 
accuracy  of  tha  maaaurad  raoaaalon  rata. 

infarrad  haating  rata  valuaa  wara  datarmlnad  for  tha  baaallna 
taat  (Run  201).  Thaaa  haating  rataa  ara  ahown  In  Flgura  3-6.  Tha  un- 
oartalnty  band  ahown  for  tha  Infarrad  haating  la  dua  primarily  to  un- 
cartalnty  In  tha  maaauramant  of  normal  aurfaoa  raoaaalon  rata.  Tha 
unoartalnty  In  normal  raoaaalon  rata  la.  In  turn,  dua  to  unoartalnty  In 
both  aurfaoa  angla  maaauramant  and  raoaaalon  rata  oomponant  (althar  axial 
or  radial) . Thla  Infarrad  haating  dlatrlbutlon  la  oomparad  with  oalorlm- 
atar  data  and  haating  pradlotlona  In  Sactlon  3.3. 


atMMt  Wk«TU  (•>  •>  W. 

Figure  3-6.  Inferred  Dlatrlbutlon  of  Nonblown  Meat  Transfer  Coefficient 
for  the  Final  Shape  of  Run  201 


3.2  SCALLOP  CALORIMETER  TESTS  (SERIES  H) 

A matal  calorimeter  replica  of  the  final  ablated  shape  of  the 
baseline  test  (Run  201)  was  oonatruoted  and  tasted  In  tha  Naval  Ordnance 
Laboratory  (NOL)  No.  8 wind  tunnel  on  February  26,  1974.  The  objectives 
of  thaaa  testa  ware  toi 

1.  Generate  heat  transfer  data  for  a calorimeter  replloa  of  an 
ablated  LTA  scallop  roughened  surface  for  comparison  with  heat 
transfer  data  Inferred  from  tha  racaaalon  rate  of  the  LTA  model. 
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2.  0«n«r«te  h*«t  transfar  data  for  a acallop  calorlmatar  modal 

ovar  a ranga  of  Raynolda  numbara  for  oomparlion  with  pradlotion 
taahnlquaa,  and  othar  rough  wall  haatlng  data. 

Xt  waa  fait  that  if  tha  calorlmatar  raaulta  eomparad  favorably  with  tha 
Infarrad  haat  tranafar  raaulta,  than  othar  camphor  noaatlp  ablation  data 
oould  bo  uaad  to  ganerata  additional  aoallop  affacta  data  for  a wide 
ranga  of  taat  oondltlona  and  noaatlp  ahapoa. 

Tha  oalorlmotor  modal  and  taat  approach  ara  daaorlbad  briefly  In 
Saotlon  3.2.1  and  tha  teat  raaulta  ara  oovarad  In  Saotlon  3.2.2. 

3.2.1  Calorlmatar  Modal  and  Teat  Approach 

Tha  purpoaa  of  tha  calorlmatar  modal  la  to  maaaura  oonvaotiva  haat 
flux  to  a aurfaoa  which  la  a rapllca  of  an  ablated  LTA  aurfaoa.  Tha 
datormlnatlon  of  haat  flux  la  baaed  on  tha  maaaurament  of  tha  tamparature 
hlatory  of  an  laolatad  tharmal  capacitance  olamant.  Haat  flux  la  datar- 
mlnod  from  thla  temperature  hlatory  by  application  of  tha  law  of  oonaar> 
vatlon  of  energy. 

Tha  calorlmatar  model  la  a rapllca  of  tha  final  ablated  ahapa  of 
Run  201  (PANT,  Sarlaa  D) . A complete  daacrlptlon  of  the  calorlmatar  la 
given  In  Volume  xill  of  thla  Interim  Report.  Tha  configuration  of 
tha  calorlmatar  modal  la  ahown  In  Figure  3>7.  The  tharmal  oapacltanea 


k ALUMINUM 


Figure  3-7.  Calorimeter  Modal  Configuration 


•laments  sr«  sagmanta  of  an  aluminum  shall  which  was  cast  from  tha  post 
test  rsplloa  of  lun  201.  Tha  sagsMnts  arc  supported  and  insulated  from 
one  another  by  a high  tamparatura  glass  filled  epoxy.  For  maasuramant 
of  tamparatura  history  tha  aluminum  sagmants  are  instrumantad  with  a 
total  of  IB  tharmooouplas.  Figure  3*8  indicates  tha  thermocouple  looa- 
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Figure  3-8.  Definition  of  Thermocouple  Locations 


During  tasting,  data  from  each  thannocoupla  ware  sampled  at  leant 
•vary  0.112  second.  Thasa  data  ware  than  curve  fit  using  a splinad 
quadratic  least  sguaras  taohniqua  to  define  tha  darivativas  of  tampara- 
tura with  tima  (dT/d6) . Whan  conduction  lossas  from  tha  isolated  thermo 
capacitance  alsmants  are  small  tha  heat  flux  can  ba  evaluated  directly 
from  tha  tamparatura  derivative  using  tha  approximation 


whars 


• . £CV  dr 

^conv  A ae 


■ convective  haat  flux,  Btu/ft*sac 

*conv 

p - material  density,  Ibm/ft* 

C ■ matarial  spacifio  haat,  Btu/lbm*R 
V ■ segment  volume,  ft* 


A ■ Mgn«nt  projsotad  araa,  ft* 
dT/d6  ■ darlvativ*  of  tamparatura  with  tlma,  *R/aao 

Analysla  prasantad  in  Voluna  xxil  Indioata  that  Equation  (3-1)  It  tuffioiantly 
aoourata  (arror  ia  lata  than  10  paroant)  if  appliad  to  timaa  laaa  than 
liO  aaoond. 


3. 2. a Taat  Kaattlta 

A aariaa  of  aavan  taata  datignatad  Taat  Sariat  H wara  parformad 
with  tha  aapliea  ealoriaatar.  Tha  taat  oonditiona  art  givan  in  Tabla  3-3. 

TABU  3-3 

TI8T  C0N01TI0M8  FOR  PART  SERIES  H HIND 
TUNNEL  TESTS  - FEBRUARY  1974* 


Run** 

No. 

Air  Supply 

Free 

Stream 

Expoaura 

Tima 

(aao) 

(paia) 

^o 

CF) 

B 

Re, 

(10*/ft) 

IS! 

65 

392 

4.95 

2.61 

HB 

370 

384 

4.99 

10.93 

B9i 

170 

375 

4.97 

5.16 

S04 

258 

379 

4.98 

7.74 

505 

369 

345 

4.99 

11.76 

13.7 

506 

675 

352 

5.00 

21.15 

12.1 

507 

503 

350 

5.00 

15.78 

12.0 

*A  tingle  oalocimatar  waa  teatad  which  waa  a replica  of 
PANT,  Sariea  D,  low  tamparatura  ablator,  Run  201,  final 

ahape . 

^All  runt  at  zero  angle  of  attack. 


During  each  taat,  ahadowgraph  photographa  wara  taken  to  define  tha  bow 
ahock  ahapa  and  datermine  whether  flew  anomaliaa  wara  occurring.  One 
of  two  typaa  of  photographa  were  taken,  either  a cloae-up  or  an  overview, 
a aanpla  of  each  typa  ia  ahown  in  Figure  3-9. 

Run  502  rapreaenta  the  moat  comparable  teat  condltiona  between 
tha  calorimeter  and  tha  LTA  taat  from  which  tha  calorimatar  waa  conatructed. 
Haaaured  early  time  (0  - 1.0  aacond)  heat  flux  lavala  for  Run  502  are 
ahown  aa  a function  of  atraam  length  in  Figure  3-10.  Certain  tharmo- 
couplaa  had  conaiatantly  auapioioua  raaponaa  and  theaa  tharmocouplaa 
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b.  Close-up,  Run  502 


Figure  3-9.  Sample  Shadowgraph  Photographs 
3-15 


Figure  3-10.  Meeiured  Heat  Flux  Levels  for  Run  502 
(Re.  - 10.93  X lOVft) 


ere  shown  with  flagged  symbols.  Also,  the  definition  of  the  projected 
area  used  in  Equation  (3-1)  for  the  stagnation  region  segment  is  uncertain 
and,  therefore,  the  heat  flux  is  also  unoertain. 

For  purposes  of  ocaqparison  with  ether  data  the  heat  transfer  coef- 
ficient is  more  applloable  than  the  heat  flux.  The  heat  transfer  coef- 
ficient was  oaleulated  frcei 


h 


^conv 

T'""Z"T 
*r  *w 


where 

h - convective  heat  transfer  coefficient 


■ boundary  layer  edge  reoovery  temperature 
T^  - wall  temperature 


(3-2) 


The  reoovery  temperature  (T^)  was  oaloulated  using  the  environmental 
routines  in  the  PANT  developed  shape  change  code  (Reference  3-4) . Com- 
parisons of  measured  heat  transfer  coefficient  with  other  data  and  pre- 
dictions are  given  in  Section  3.3. 
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3.3  DATA  COMPARISONS 


Thr««  Mts  Of  rolovant  rough  wall  haat  transfer  data  are  compar- 
able i 

• Heat  transfer  rates  inferred  from  the  measured  recession  rate 
for  Run  201 

• Heat  transfer  data  measured  from  the  scallop  replica  calorim- 
eter model  (Series  H) 

e Heat  transfer  data  measured  for  sandgraln  roughened  calorim- 
eter models  (Ssrles  A) 

A comparison  of  the  inferred  heat  transfer  rates  and  scallop  replica  data 
la  given  In  Section  3.3.1.  Comparisons  of  the  scallop  calorlmetar  data 
with  the  PANT  correlation  developed  from  the  Series  A sandgraln  roughness 
calorimeter  data  are  discussed  In  Section  3.3.2.  Conclusions  regarding 
these  data  comparisons  are  presented  In  Section  3.3.3. 

3.3.1  Comparison  of  inferred  Heat  Transfer  Rates  and  Scallop  neplloa 


Since  the  scallop  replica  model  was  constructed  of  materials  which 
degrade  at  higher  temperatures « all  Series  H tests  were  aonduoted  at  sup- 
ply temperatures  below  400*F.  The  comparable  LTA  model  (Run  201)  was, 
however,  tested  at  a supply  temperature  of  S85*F.  in  order  to  compare 
the  Inferred  heating  rates  from  the  ablation  tost  with  the  calorimeter 
data,  the  Inferred  heating  was  corrected  to  account  for  differences  In 
supply  temperature  and  pressure  as  described  in  Volume  XIII.  Figure  3-11 
shows  the  adjusted  Inferred  heating  distribution  and  the  comparable 
calorimeter  data  (Run  502,  Re^  ■>  10xl0*/ft) . The  boxed  areas  around  the 
calorimeter  data  Indicate  the  spread  In  the  data  on  a given  ring  and  the 
stream  length  covered  by  each  ring.  Because  the  rings  are  Insulated  from 
one  another  the  measured  heat  flux  represents  an  average  of  the  local 
heat  flux  over  the  ring.  The  data  are  in  general  agreement  with  Infer- 
red heat  coefficients  from  the  LTA  data.  The  calorimeter  data  on  the 
forecone  (Rings,  2,  3,  and  4)  agree  quite  well  with  the  Inferred  haat 
transfer  but  the  aft  cone  data  (Ring  5)  Is  higher  than  the  Inferred. 


9rB4UU  ULhatU 

Flgdr*  3-11.  Comparison  of  infsrred  and  Msasurad  Hast  Transfer  coeffloient 

Distribution  for  the  Baseline  Reynolds  Number  Condition  (Run  502) 

Indeed,  the  mean  of  the  oalorimeter  data  inoreases  relative  to  the  mean 
of  the  inferred  values  as  a function  of  streamwise  distanoe.  It  may  be 
speculated  that  this  trend  is  assooiated  with  the  accumulating  effects 
of  ablation  products  on  the  boundary  layer  heat  and  mass  transfer.  For 
the  ablating  case,  mass  addition  to  the  boundary  layer  over  the  foreeone 
would  reduce  mass  transfer  (and,  hence,  inferred  heat  transfer)  at  aft 
cone  locations.  No  calculations  wore  performed  to  assess  this  specula- 
tion. 


3.3.2  Comparison  of  Scallop  Renlica  Date  with  Sandarain  Roughness 
Effects  bate  and  Correlations 

Under  the  PANT  program  extensive  testing  of  sandgrain  type 
roughened  calorimeters  was  performed  (Secies  A tests  ~ Reference  3-5) . 

The  data  from  these  tests  wes  used  to  develop  a rough  wall  heating  corre- 
lation applicable  to  sandgrain  type  surfeoo  roughness  (Reference  3-6) . 

In  this  correlation  the  roughness  heating  augmentation  factor 


rough 

'smooth 


(3-3) 
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is  corrslatsd  varsus  s parameter  which  is  approximately  proportional  to 
the  ratio  of  the  roughness  height  to  the  laminar  sublayer  thickness 


( 


R«. 


(3-4) 


where 


k " sandgraln  roughness  height 
- laminar  sublayar  thlckneaa 

Rej^  “ ■ - boundary  layer  edge  Reynolds  number  based  on 

^e  roughness  height 


Ciig  ■■  smooth  wall  Stanton  number 

" local  boundary  layer  edge  temperature 
■ local  wall  temperature 


A comparison  of  the  Series  K calorimeter  data  with  the  PANT 
rough  wall  heating  correlation  Is  shotm  In  Figure  3-12.  The  value  of 
for  each  data  point  was  calculated  by  dividing  the  measured  heat 
transfer  coefficient  by  the  smooth  wall*  turbulent  PANT  code  prediction. 
The  error  bands  on  each  represent  the  spread  In  the  data.  The  abscissa 
is 


R«. 


(3-5) 


where 


Re. 


0 u 

^e  e 


Therefore,  the  PANT  correlation  appears  as  a series'  of  parallel  lines, 
one  for  each  roughnesj  height.  The  value  of  the  abscissa  parameter  for 
each  data  point  was  taken  from  the  smooth  wall  predictions.  Different 
symbols  are  shown  for  each  wind  tunnel  run.  The  figure  Indicates  that 
Is  roughly  1.3  on  the  forecone  ring  (Ring  2),  increases  to  1.6  on 
Rings  3 and  4,  and  further  Increases  to  approximately  2.0  on  Ring  5. 

This  observation  Is  consistent  with  the  fact  that  scallop  dimensions  are 
greater  on  the  aft  cone  than  on  the  forecone.  The  point  of  Interest  Is 
that  the  heating  increase  factor  on  each  ring  Is  not  sensitive  to  the 
test  Reynolds  number.  The  measured  heating  rates  scale  closely  with 
changes  In  smooth  wall  turbulent  heating.  This  demonstrates  that  the 
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xefwar' 


Heat  AuyKntatiow  Factor  (t^ 


•oallopi  affaot  h««t  trutsfar  in  a di££ar«nt  way  than  do  sandgrain  rough- 
naiB  alamanta  and  tanda  to  oon£lnn  tha  apaoulation  atatad  in  Ra£eranoa 
3-6  that  tha  PANT  roughnaaa  a££aota  oorralatlon  ia  not  diraotly  appllo- 
abla  to  aoBllopad  aur£aeaa.  Indaad,  ainoa  aoallopa  £onn  aa  a raault  o£ 
£low/taatarial  ablation  intaraationa»  it  would  not  ba  axpaetad  that  a 
oorralatlon  o£  aandgrain  haating  data  would  alao  oorralata  aoallop  haat- 
ing  data. 

Sinea  tha  augmantation  ia  inaanaltiva  to  Raynolda  numbar  condition, 
ona  might  oonjaetura  that  aoallopa  a££aet  haating  mainly  through  ohangaa 
in  tha  invlaoid  flow  rathar  than  through  ohangaa  in  tha  boundary  layar. 
Tha  ahadowgrapha  ahown  in  Flgura  3-9  confirm  that  aignifioant  aurfaoa/ 
ahook  layar  intaraotiona  did  ooour  in  tha  oalorimatar  taata.  Similar 
intaraotiona  wart  obaarvad  in  ahadowgrapha  from  tha  Sariaa  D and  Sariaa 
Z-LTA  taata. 

Another  intaraating  point  ia  that  tha  aoallop  aurfaoa  produoaa 
laaa  haating  augmantation  than  a aandgrain  aurfaoa  having  a aimilar 
roughnaaa  height  taated  at  tha  aama  Raynolda  numbar  oondition.  Thia 
point  ia  llluatratad  in  Figure  3-13  whara  tha  aonie  point  augmantation 
faotor  from  tha  aoallop  oalorimatar  (Ring  2)  ia  compared  to  oorraapond- 
ing  data  from  tha  Sariaa  A,  40  mil,  aandgrain-typa  roughnaaa  oalorimatar 
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Figure  3-13.  Comparison  of  tha  Soallop  Calorimeter  Data 
With  Sandgrain  Type  Calorimeter  Data 
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modal.  Tha  roughnaai  maasuramants  ihown  in  Saotlon  3.1.1  indloata  a 
nominal  roughnaaa  halght  of  about  30  mils  for  tha  aonio  ragion  (i.a., 
foraeona)  of  tha  aoallop  oalorlmatar.  Although  tha  roughnaaa  halght  is 
similar  bstwasn  tha  two  modals  (30  mils  varsus  40  mils)  tha  sandgrain 
haating  factor  is  significantly  higher  than  tha  oomparabla  scallop 
calorimatar  result.  A taohniqua  described  in  Rafaranea  3-7  to  datarmina 
equivalent  sandgrain  roughnass  dimensions  from  scallop  gsomatry  data 
was  investigated  in  an  attempt  to  rationalisa  tha  diffaraneas  shown  in 
Figura  3-13.  Although  judicious  specif ioation  of  scallop  geometry  param- 
atars  can  explain  tha  relative  values  of  haating  factors « tha  different 
trends  with  Reynolds  number  condition  could  not  be  axplalnad.  To  modal 
the  observed  trand«  tha  equivalent  sandgrain  would  have  to  be  a function 
of  test  condition. 

3.3.3  Conclusions  and  Reocamandations 


3. 3. 3.1  Conclusions 

Tha  following  conclusions  are  dratm  from  observation  of  the  data 
and  comparison  with  predictions. 

a The  heat  transfer  to  the  scalloped  surface  calorimatar  was 
significantly  greater  than  corresponding  smooth  wall  predic- 
tions . 

a Tha  measured  heat  transfer  rates  agree  generally  with  heating 
rates  inferred  from  the  corresponding  LTA  model;  differsnoes 
suggest  possible  affects  of  ablation  on  the  inferred  heat 
rates . 

a The  roughness  augmentation  factors  (K^)  for  each  scallop 

calorimeter  ring  do  not  scale  with  the  PANT  sandgrain  rough- 
ness effects  correlation  parameter 


The  roughness  augmentation  factor  (X^)  distribution  around 
the  body  is  nearly  independent  of  free  stream  Reynolds  number, 
but  increases  with  increasing  scallop  roughness  height. 

The  scallop  surface  on  the  Series  H calorimeter  produces  less 
heating  augmentation  than  a sandgrain  surface  having  a similar 
physical  roughess  height. 


If 


• The  •qulval«nt-*andgrain*rott9hn«sa  fomulatlon  of  Rafsranoa 
3*7  la  not  adaquata  to  axplain  diffaranoaa  batwaan  aandgrain 
oalorlmatar  and  acallop  oalorinatar  data. 

3. 3. 3. 2 Raoonnandations 

Tha  following  raaonnondatlons  raault  from  tha  aoallop  atudiaa 
parformad  undar  thla  program 

a Furthar  aasassmant  of  tha  affaets  of  upatraam  ablation  on 
Infarrad  haatlng  rataa  ahould  ba  parformad  to  eatabllah  con- 
olualvaly  that  LTA  ablation  data  oan  ba  uaad  to  Infar  noaa- 
tip  flow  phanomana. 

a Infarrad  haatlng  rataa  ahould  ba  darlvad  from  a largo  aot  of 
LTA  taat  raaulta  to  provlda  data  for  acallop  affaota  atudiaa. 

a Corralatlona  of  acallop  affaota  on  haat  tranafar  ahould  be 

darlvad  aa  funotlona  of  Invlaeld  and/or  vlacoua  flow  proportlaa. 

a Tha  validity  of  tha  corralatlona  darlvad  from  LTA  wind  tunnol 
data  ahould  be  verified  through  oomparlaona  to  high  preaaure 
ablation  data  of  actual  noaetlp  matarlala 
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8ECTZ0K  4 

SERIES  Z WZND  TUNNBI.  TESTS 

Th«  SariM  Z Low-Tamparatura  Ablator  (bTA)  wind  tunnal  taata  wara 
oonduotad  In  Tunnal  Mo.  8 at  tha  Naval  Ordnanoa  Laboratory  (MOL)  from 
March  27  to  April  4«  1974.  Tha  ganaral  purposa  of  thaaa  taiti  wai  to 
ganarata  LTA  shapa-ohanga  data  to  datarmlna  tha  affaot  of  ralavant  nosa- 
tlp  paramatara  on  tha  formation  and  axtant  of  Irragular  shapaa.  A sum- 
mary of  apaoifio  taat  objaotlvaa  and  approach  ia  glvan  in  Saction  4.1i 
aamplaa  of  tha  raaulting  data  ara  givan  in  Saction  4.2)  and  data  tranda 
ara  aummarlaad  in  Saction  4.3.  A datailad  praaantation  of  tha  taat  data 
are  alao  givan  in  Volume  XZV  of  thla  report  aariaa. 

4.1  TEST  OBJECTIVES  AMD  APPROACH 

The  objactlvaa  of  thla  axparimantal  program  wara  aa  followai 
a Evaluate  tha  affacta  of 


- Raynolda  number 

- Tranaition  control 
— Surface  roughnaaa 

- Model  alsa 

- Initial  geometry 

- Blowing  rata  and  wall  temperature  ratio 

on  noaatip  aharpening  and  the  formation  and  axtant  of  Irragular 

ahapaa. 

a Ganarata  accurate  racaaaion  data  to  evaluate  heat  tranafer 
rates  to  scalloped  aurfacaa. 

a Ganarata  shape-change  and  shock-ahapa  data  to  evaluate  the 
phenomena  associated  with  the  growth  and  degradation  of  small 
nose  protubaraneaa . 

a Generate  high-frequency  acceleration  data  to  evaluate  tha 
effects  of  oscillating  shocks. 

a Simulate  flight  shape-change  raaponse  in  a wind  tunnal. 


■ V 


I- 


ThM«  objtotlvai  vara  addraaaad  by  axposing  30  oaaphor  lta  modals  of 
0 diffarant  initial  gaonatriaa  to  various  fraa  straam  conditions  in  tha 
MOL  No.  8 hyparsonio  wind  tunnal  at  a Naoh  nusibar  of  5.  A typical  modal 
gaonatry  is  shown  in  Vigura  4-li  and  tha  tast  matrix  is  givan  by  Tabla 
4*1.  Tha  larga  amount  of  data  acquired  to  acconplish  tha  tast  objactiva 
inoludad  18  nm  and  35  sat  shadowgrapha  # 18  am  movias^  35  an  and  70  nn 
photographs!  and  high-fraquanoy  aooalaration  aaasuramants . 

4.2  8ANPII  DATA  NB8ULT8 

Tha  data  of  primary  intaraat  for  this  tast  sarias  ara  tha  shapa* 
ohanga  historian  of  tha  LTA  modals.  This  Infomatlon  Is  availabla  in 
raw  photographio  data  consisting  of  18  an  aldaviaw  shadowgraph  movies ! 

35  am  sidaviaw  shadowgraph  pulaad  photographs#  35  mm  overhead  pulsed 
photographs#  and  70  an  underside  pulsed  photographs,  To  raduca  the 
raw  data  to  an  easily  usable  form#  tracings  ware  mada  from  tha  35  mm 
film  to  record  nodal  shape  profile  historian.  An  axampla  of  this#  for 
Run  805#  is  shown  in  Figure  4-2.  In  addition#  further  reduced  data 
related  to  shape-change  history  ara  typified  by  tha  stagnation  point 
axial  raoassion  history  given  for  Run  805  in  Figure  4-3. 

Data  of  secondary  intaraat  for  Sarias  I inoludad  boundary  layer 
transition  locations  and  subaaquant  sharpening  phanomana.  Initial  transi- 
tion locations  ware  determined  from  tracings  of  35  iwn  overhead  and  70  mm 
underside  photographs#  as  shown  in  Figure  4-4#  and  ware  subsequently  re- 
duced to  (S#4)  coordinates  as  shown  by  Figure  4-5.  Both  qualitative  and 
quantitative  transition  data  are  also  available  from  tha  70  mm  front- 
view  photographs  and  16  mm  overhead  high-speed  movies#  as  demonstrated 
by  Figures  4-6  and  4-7#  respectively.  Information  relative  to  sharpening 
can  be  taken  from  most  of  the  film  sources#  but  is  given  most  clearly 
by  the  70  mm  underside  photographs#  of  which  examples  are  presented  in 
Figure  4-8. 

Additional  data  of  Interest  inoludes  those  related  to  unsteady 
flow  fields  with  oscillating  shocks  on  Irregular  shapes.  This  area  of 
interest  was  addressed  primarily  through  the  usage  of  a vary  high  spaed 
(20 #000  frames/second)  Hycam  16  mm  camera  taking  shadowgraph  movies#  and 
a three-axis  high-frequency  accelerometer  system  mounted  on  the  sting. 

The  Hycam  was  triggered  manually  during  the  tost  when  accelerometer 
output  increased.  Figure  4-9  demonstrates  typical  accelerometer  re- 
sponses during  periods  of  unsteady  flow.  Oscillating  shocks  during 
those  periods  are  also  visible  on  the  16  mm  and  35  mm  shadowgraph  film; 
an  example  of  the  latter  is  shown  in  Figure  4-10. 
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TABLE  4-1 

TEST  CONDITIONS  FOR  PANT  SERIES  I TESTS 
MARCH/APRIL  1974 


Titt  eondltloM 


Prtt  StriMi 


fiiplior  Ihdtl 


1.8*  8(1  Sphart-coM 

1.8*  R||  SpharfCOM 

1.3*  Up  Splitraiiema 

1.8*  Up  Spliar^coMt  partlcla-lMdad 

1.8*  8p  SpHara^oM 

0.78*  lp\Spliart>ceM 

1.8*  Rp  SpiMfa-eont  with  groove 

1.8*  Rp  Sphoro^ono.  partlelo-loidod 

0.78*  Rp  Sphort-eono 

1.8*  Rs/80°/8^  Rieenic 

1.8*  Rp  Sphoro-cono.  ptrtlelo-loadod 

0.78*  Rp  Sphtre*eono 

1.8*  Rs/TS^/O**  Riconie 

1 .8*  Rp  $ptioro<OM  with  3 groovtt 

1.8*  Rg  LMlMr>h1unt 

1.8*  Rs/60®/8®  Bicoiilc 

1.8*  Rp  Sphoro-com 

1.8*  Rp  Sphoro'cono 

0.75*  Rp  Sphtrt>cont 

3.8*  Rj  LwlMr>b1unt 
2.5*  Rp  Sphtrt-cem 
1.5*  Rg  Lawinor-hlunt 
1.8*  R5/75V8*  Bleonic 
1.5*  Rp  Sphtro>cono 
2.5*  Rp  Sphert-cono 
1.5*  1^  LMrtntr«b1unt 
2.8*  Rg/S5°/6°  BIconIc 
0.78*  Rp  Sphoro-cono 
0.75*  Rp  Sphort-eoM 
3.8*  Rj  Lw1nar>b1unt 
3.8*  R«  Loiilmr-blifflt 
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(SM) 


8B.63 

10B.47 

172.83/ 

308.43 

63.49 

118.74 

172.91 

123.15 

126.28 

328.48 

253.84 

20.84/ 

69.04 

122.83 

13.47 


130.48 

242.75 

67.73 

192.33 

139.04 

91.11 

116.03 

153.52 

231.08 


Dual  Ro.  test  - valMi  glvoo  art  for  first  portlon/soeont  portion 
Approxlnato  tliaa  onlyt  fllw  data  Ineowpltta 


Figure  4-6.  Typical  70nn  frontview  photographs  (Run  805;  1.5-inch 
Rjj  sphere-cone;  Re^  = 6.56  x 10*/ft,  = 987“F). 


Model  exposure  time 
approximately  O.S  second 


Figure  4-7. 


Example  of  Transitional  Gouge  Development  from  16mm  Overhead 
Camera.  (Run  805;  1.5-1nch  Sphere-Cone;  Re*  ■ 6.56  x 10*/ft, 
Tq  ■ 987  F) 
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Figure  4-9.  Example  of  Acceli 
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4.3  PRELIMINARY  DATA  TRENDS 

Tha  primary  objectlva  of  thla  tact  aarlai  ai  dasorlbad  In  Saetlon 
4.1,  was  to  avaluata  the  affects  of  six  significant  parameters  on  tha 
formation  and  extent  of  irregular  shapes.  Preliminary  examination  of 
the  data  leads  to  ths  following  observations  in  regard  to  each  parameters 

e Reynolds  number  - As  ehotm  in  Figure  4-11,  free  stream  Reynolds 
number  has  a strong  effect  on  shape  change.. 

e Transition  control  - Boundary  layer  trip  grooves  used  on  Ru'ia 
807  and  814  did  not  produoe  completely  symmetric  transition 
at  the  low  Reynolds  number  conditions  tested.  However,  shape 
profile  development  in  these  runs  was  different  than  exhibited 
by  nongroved  models. 

e Surface  roughneea  - Glass  particles  imbedded  in  the  models 

for  Runs  804,  808,  and  811  had  a strong  affect  on  shape  change, 
apparently  impeding  the  formation  of  the  irregular  shape  regime. 

e Model  else  - Figure  4-12  illustrates  that  decreasing  the  modal 
else  (nose  radius)  causes  an  increase  in  the  Reynolds  number 
required  to  attain  a given  shape  regime. 

e Initial  geometry  - A comparison  of  post-sharpening  and  final 
shapes  for  models  with  various  initial  geometrlss  is  shown  in 
Figure  4-13.  It  appears  that  there  are  only  slight  dif- 
ferences in  final  shape  among  the  different  initial  geometry 
models.  However,  different  shapes  led  to  different  times  for 
sharpening. 

e Blowing  rate  (B‘)  and  wall  temperature  ratio  (Tw/T»)  - Tests 
run  at  low  supply  temperatures  (T^  ■ 500*F  compared  to  the 
baseline  T^  - 1000*F)  result  in  lower  B'  and  higher 
The  combined  effect,  as  shown  in  Figure  4-14,  is  an  Increase 
in  the  Reynolds  number  required  to  attain  a given  shape  regime. 

An  additional  significant  trend  was  noted  in  regard  to  unsteady 
flow  fields;  namely,  unateady  flow  phenomena  existed  for  portions  of 
about  two-thirds  of  the  runs,  and  for  portions  of  virtually  all  runs 
which  developed  irregular  shapes.  Shape  change  after  flow  Instabilities 
started  tended  to  produce  shapes  which  damped  the  instabilities.  Cycles 
of  steady  to  unsteady  flow  were  observed. 
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Figure  I'll.  Raynolda  Number  Effect  on  Model  Plnel  Shape 

(All  1.5-Znoh  Sphere-Conee » Nominal  Tg  ■ lOOO'F) 
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(b)  Nominal  ■ 5 k 10*/ft 

Flgura  4-13.  Comparison  of  Poat-Sharpaning  and  Final  Shapaa  for  Diffarant 
Initial  Gaomatry  Modala 
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SERIES  J WIND  TUMMBL  TESTS 

During  th«  PANT  program  two  ealorlmatar  taat  aariaa  wara  parformad 
to  aatabllah  oonditiona  raqulrad  to  promota  boundary  layar  tranaltlon  on 
rough  axlaymmatrlo  blunt  bodiaa.  Tha  first  taat  program  la  rafarrad  to 
aa  Sarias  A;  tha  aaoond  as  Sarlaa  J.  Tha  raaults  of  Sarins  A ara  doeu- 
msntad  in  Rafaranoaa  5-1,  5-2,  and  5-3  which  daaoriba  tha  taat  data, 
roughnaaa  augnantad  haating  analysis, -and  boundary  layar  transition  analy- 
sis, raspaotivaly.  In  this  ssotlon  a summary  of  tha  Sarias  J tast  pro- 
gram is  givan.  This  inoludas  a disousaion  of  tha  sxparimantal  tast  ob- 
jaotlvas  and  approach  (Section  5.1),  aa  >wall  as  a briaf  ovarviaw  of  tha 
significant  tast  rasulta  (Saction  5.2). 

5.1  BXPERZNBNT  OBJECTIVES  AND  APPROACH 

Tha  taat  objactivas  wara  dafinad  baadd  on  tha  data  dafioianoias 
which  wara  notad  at  tha  complntion  of  tha  Sarias  A oorralation  efforts. 

Tha  primary  tast  objactivas  warat 

a To  avaluata  tha  affaot  of  body  siss  on  transition  onsat  and 
location 

a To  avaluata  tha  affact  of  body  shape  on  transition  onset  and 
location 

Sacondary  objactivas  wara  to  obtain  additional  data  on  tha  affaot  of  tha 
wall  tamparatura  ratios,  * on  rough  wall  boundary  layar  transition 

and  to  obtain  additional  data  ralativa  to  rough  wall  heating  augmentation. 

As  can  be  inferred  from  tha  taat  objectives  tha  primary  dafioianoy 
of  the  Sarias  A program  was  tha  lac)c  of  transition  data  on  different  sisa 
and  shape  bodies.  Although  the  relative  importance  of  these  two  param- 
eters was  not  known,  the  Sarias  A data  did  Indioata  tha  possibility  of 
those  paramatars  influoneing  boundary  layar  transition.  This  conclusion 
was  primarily  based  on  tha  axparlmantal  observation  that  transition  on 
tha  roughanod  calorimeter  models  never  ooourrad  batwasn  tha  sonic  point 
and  tanganey  point.  It  was  postulated  that  in  this  ragion  trusitlon  is 
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■uppraasad  by  tha  indapandant  or  ooabinad  affaota  of  tha  nomal  praaaura 
gradlant  asaoolatad  with  eantrlfugal  aeealaratlon  and/or  tha  atraamwiaa 
praaaura  gradlant.  For  aphara  oona  modala,  both  praaaura  gradlanta  ara 
proportional  to  tha  Invaraa  of  tha  noaa  radiua.  Manoa,  body  alia,  in 
partloular  noaa  radiua,  was  a paramatar  which  daaarvad  further  conaldara- 
tlon. 

Tha  Sariaa  K taat  program  only  partially  addraaaad  tha  Inqportanoa 
of  body  ahapa  on  boundary  layar  tranaltlon.  Dlffaranoaa  in  body  ahapa 
raflaot  dlffaranoaa  in  tha  praaaura  gradlanta  and,  hanoa,  oould  ba  ax- 
paotad  to  influanoa  boundary  layar  tranaltlon.  Furtharmora,  body  ahapa 
oould  influanoa  tranaltlon  through  tha  affaota  of  boundary  layar  apraadlng. 

Six  modala  wara  aalaotad  to  aooompllah  tha  above  atated  objaotlvaa. 
Tha  modal  oonflguratlona  ara  illuatratad  In  Flguraa  5-1  and  5-2.  Four 
of  tha  modala  wara  aphara^onea  having  noaa  radii  of  0.75,  1.5,  2.5  , and 
3.5  Inohaa.  Tha  half  oona  angle  for  thaaa  modala  wan  8*.  The  two  ramaln- 
Ing  models  wara  a 60*  bloonio  modal  and  a laminar  ablated  shape  modal. 

Tha  oontinuous  thick  akin  ealorlmatry  taohnlqua,  used  suooass- 
fully  during  tha  Sariaa  A teat  program,  waa  aalaotad  for  obtaining  haat 
transfer  measuramants.  The  Series  J models  wara  machined  from  Niokal  200 
to  a nominal  wall  thioknass  of  0.080  Inch.  Tha  oalorlmater  shells  were 
inatrumantad  with  36  gauge  (5  mil)  Chromel-Alumal  thermocouples.  The  num- 
ber of  tharmooouplas  dapandad  upon  tha  alia  of  tha  modal)  tha  larger 
modala,  such  aa  tha  3. 5- Inch  sphara-oona  modal,  had  78  tharmooouplas 
whereas  tha  smallest  modal  had  39  tharmooouplas.  All  of  tha  models 
except  for  the  Series  A 2.5-inoh  nose  radius  model  were  grit  blasted  to 
obtain  a uniformly  distrlbutsd  paak-to-vallay  roughness  of  3.5  mils.  Tha 
Series  A modal  had  a roughness  of  3 mils.  Additional  details  relative  to 
the  models  may  be  found  in  Volume  XV  of  this  report. 

Tha  Series  J test  program  was  conducted  in  the  Naval  Ordnance 
Laboratory  (NOL)  Tunnel  No.  8 from  February  22,  1974  to  February  26,  1974. 
This  faoility  was  ohoaan  aarly  in  tha  PANT  program  baaed  on  its  capability 
of  providing  hyparsonio  flow  and  high  sonlo  point  unit  Reynolds  numbers 
(Rafarenoa  5-4) . Tha  free  stream  Naoh  number  for  tha  Series  J tests 
waa  5. 


*Tha  2.5-inoh  modal  waa  tasted  during  Sariaa  A and  was  reused  for  Series 
J. 
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EXPERIMENTAL  RESULTS 


The  Seriee  J teit  program  suooaesfully  accomplished  all  stated 
objectives.  The  effects  of  body  else  (nose  radius)  and  oody  shape  on 
boundary  layer  transition  wore  adequately  characterized.  A total  of  40 
calorimeter  model  entries  were  made  In  the  NOL  wind  tunnel,  spanning  a 
sonic  point  unit  Reynolds  number  range  from  0.42  x 10*/ft  to  7.4  x 10*/ft. 
A summary  of  the  Series  J test  conditions  is  given  in  Table  5-1.  In 
this  section,  a brief  aumnary  Is  given  of  the  data  reduction  procedures 
as  wall  as  some  of  the  mere  salient  results. 

The  Series  J data  were  reduced  using  the  techniques  developed 
during  the  Series  A test  program.  These  techniques  include  temporal 
curve  fitting  of  the  individual  thermocouple  ten^erature  data,  calcula- 
tion of  the  local  incident  heat  flux  accounting  for  straamwlse  conduc- 
tion, and  calculation  of  the  local  heat  transfer  coefficient  based  on 
recovery  temperature.  These  techniques  are  documented  In  Referenes  5-3 
and  Volume  XV  of  this  report. 

Using  these  data  reduction  techniques,  heat  transfer  coefficients 
can  be  presented  either  as  a function  of  meridian  location  or  as  a func- 
tion of  time.  For  purposes  of  this  report,  heat  transfer  coefficient  dis- 
tributions presented  as  a function  of  meridian  location  will  be  referred 
to  as  early  time  data.  This  reflects  the  fact  that  the  heat  transfer 
coefficients  are  evaluated  at  one  time,  typically  1 second.  Similarly, 
heat  transfer  coefficient  distributions  presented  as  a function  of  time 
will  be  referred  to  as  late  time  data.  In  this  case,  the  heat  transfer 
coefficients  are  calculated  at  various  times  which  correspond  to  varlour 
wall  temperature  ratios  during  the  test.  The  handling  of  streamwlse  con- 
duction Is  slightly  different  for  early  time  as  opposed  to  late  time 
(refer  to  Reference  5-3  or  Volume  XV).  As  Indicated  In  the  Appendix 
of  this  rsport,  late  time  data  Is  used  In  correlating  the  transition  data. 
This  Is  necessary  since  rough  wall  boundary  layer  transition  Is  a strong 
function  of  the  wall  temperature  ratio. 

Figure  5-3  Illustrates  typical  early  time  heat  transfer  coefficient 
(HTC)  values  as  a function  of  nondimens ional  streamwise  length  (S/R^])  for 
the  1.5- inch  nose  radius  sphere  cone  model,  the  blconlc  model,  and  the 
laminar  ablated  shape  model.  Several  Interesting  trends  can  be  noted  from 
this  flgurs.  First,  the  lach  of  data  scatter  at  any  streamwlse  location 
Indicates  the  consistency  and  accuracy  of  the  calorimeter  data.  Second, 
all  of  those  runs  to  a certain  extent  show  the  existence  of  boundary  layer 
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TABLE  5-1 

ACTUAL  TEST  CONDITIONS  FOR  PANT  SERIES  J WIND  TUNNEL  TESTS 


Run 

IMn 

M. 

No. 

II 


No«t1  MKrinton^ 


Z.i  In.  R|,  Spntra/eom 


3.S  In.  R||  SpUnri/eonn 


.7S  In.  R||  Sphnrn/eont 


1.S  In.  Du  Sphnn/eniM 


llconict  Rj  ■ 1.S  In. 


LwlMr  Militnd  Snapt. 
Rj  > 3.1  In. 


All  runt  at  itro  ang1a*ef-attack 
%11  noAtl  grit  blattad 


90*0  ! oo*«P 


(a)  1.5-lnch  Sphere-Cone  Model 

Figure  5-3.  Typical  Early  Time  Heat  Transfer 
Coefficient  Distributions 
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ralaminariiAtlon  in  tha  oornar  ration  of  tha  nodals.  Thia  phanonana 
was  also  notad  in  tha  Sarias  h tasts^  and  has  also  baan  obaarvad  on  gra- 
phitic nodals  taatad  in  tha  SO  MW  AFFDL  Are  Facility  (Rafaranoa  5-5) . 

Third,  tha  transitional  haat  transfar  ooaffioiant  distributions  ara 
similar  for  tha  sphara  oona  and  laminar  ablatad  ahapa  modal  and  quita 
diffarant  for  tha  sinpla  bioonic  nodal. 

Pigura  5-4  oomparas  turbulant  haat  transfar  ooaffioiant  distribu- 
tions for  thraa  of  tha  four  sphara  oona  modals  tasted  during  Sarias  J. 

Tha  nominal  sonio  point  unit  Raynolda  numbar  for  thaaa  runs  is  1.75  x 
10*/ft.  This  figura  illuatratas  that  boundary  layar  transition  is 
affaotad  by  tha  sisa  of  tha  modal}  howavar,  tha  affact  is  not  larga. 

Indaad,  as  is  indioatad  in  Saotion  6.1,  tha  Sarias  J transition  data 
vfara  found  to  oorralata  raasonably  wall  using  praviously  darivad  sarias 
A transition  oritarion.  This  critarion  indioatas  a nose  radius  dapandanoy 
of 

Figura  5-5  illuatratas  haat  transfar  ooaffioiant  distributions 
for  tha  bioonio  modal  as  a function  of  Raynolda  numbar.  Tha  most  intar- 
astlng  obsarvatlon  is  tha  Inoraaslng  trand  of  tha  haat  transfer  coef- 
fiolant  along  tha  fora  oona  of  tha  modal.  This  trand  was  also  observed 
in  tha  Series  A teat  program  on  tha  simple  bioonio  modal  (Rafaranoa  5-6). 
Boundary  layar  solutions  performed  for  both  of  these  shapes  indicate 
that  tha  haat  transfar  on  tha  fora  oona  is  very  sensitive  to  tha  bound- 
ary layer  edge  entropy. 

Some  examples  of  lata  time  data  raduotion  ara  shown  in  Figura  5-6 
for  the  3.5-inch  sphara  oona  model,  the  bioonio  model,  and  tha  laminar 
ablatad  shape  modal.  This  figure  presents  haat  transfar  coefficient 
distributions  at  various  times  during  tha  test.  The  models  heat  up  dur- 
ing the  test,  and  different  times  correspond  to  different  wall  temperature 
conditions.  That  is,  wall  temperature  ratio  increases  as  the  time  increases. 

Laminar  haat  transfer  coefficient  distributions  are  shown  as  the 
solid  line  in  Figure  5-6.  These  predictions  wars  made  using  the  PANT 
shape  change  code  including  the  modification  described  in  Section  6.2. 

The  solutions  were  performed  at  a wall  tamparature  corresponding  to  1 
second  tast  tima;  however,  laminar  haat  transfar  is  Insensitive  to  the 
wall  tamparature  ratio  and,  hence,  these  represent  laminar  heat  transfer 
over  tha  Indicated  test  times. 

comparison  of  tha  data  with  tha  theoretical  laminar  distribution 
illustratas  the  strong  dependency  of  boundary  layar  transition  on  wall 
temperatura  ratio.  For  all  tha  models  shown  in  Figura  5-6,  the  transition 
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Figure  5-4.  Comparison  of  Sphere-Cone  Turbulent  Heat  Transfer 
Coefficient  Distribution  Data 
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(a)  3.5-lnch  Sphere  Cone  Model 

Figure  5-6.  Typical  Late  Time  Heat  Transfer  Coefficient  Diatributione 
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(c)  Laminar  Ablatad  Shapa  Modal 
Flgura  5-6.  Concludad 


5-15 


1 


location  olaarly  movaa  downatraam  as  tha  tast  tina  Inoraaaaa  or,  mora 
approprlataly,  aa  the  wall  tanparatura  ratio  inoraaaaa.  Thia  raault 
waa  adaquataly  invaatigatad  in  tha  tariaa  A taat  program  on  sphara-oona 
nodala  and  tha  Sariaa  J taat  raaulta  oonfim  tha  affaot  ot  tha  wall  tam> 
paratura  ratio  for  tha  aphara-oona  modala  aa  wall  aa  for  tha  bioonio  and 
laminar  ablatad  ahapa  modal. 

In  aumnary,  tha  Sariaa  J taat  program  haa  providad  a aubatantial 
data  baaa  from  which  tha  af facta  of  body  aiaa  and  ahapa  on  rough  wall 
boundary  layar  tranaition  can  ba  dataminad. 
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SECTION  6 


IMPROVED  N08BTIP  PERFORMANCE  ANALYSIS  TECHNIQUES 

At  th«  complation  of  testing  under  PANT  Teak  4.2.5,  aaverel 
anelyele  teaks  were  performed  to  upgrade  the  oomputer  code  developed 
under  previous  PANT  tasks  (4.2.1  through  4.2.4).  The  efforts  consisted 
of  either  the  correlation  of  PANT  test  data  to  determine  shape  change 
response  regimes  or  the  evaluation  of  improved  computational  tools  for 
nosetip  design  analysis.  Response  trends  from  the  respective  test  series 
are  presented  in  the  test  summary  sections.  In  this  section,  improved 
modeling  procedures  are  presented  and  results,  using  the  improvements, 
are  compared  to  results  from  previous  calculations  for  seven  baseline 
flight  environment  oases.  Section  6.1  sunmariies  recommended  modifica- 
tions to  the  PANT  rough  wall  transition  criterion.  Section  6.2  describes 
a modified  technique  for  computing  turbulent  heat  transfer  on  a reentry 
vehicle  nosetip.  Section  6.3  discusses  numerical  techniques  used  to 
calculate  nosetip  shape  and  shape  change  response.  The  comparative 
calculation  results  are  presented  in  Section  6.4. 

6.1  UPDATED  TRANSITION  CRITERIA 

A detailed  discussion  of  the  PANT  rough  wall  transition  studies 
is  presented  in  the  appendix  to  this  report.  The  primary  results  are 
outlined  in  the  following  summary. 

Two  experimental  programs  (Series  A and  J)  carried  out  at  • 5 
in  NOL  wind  tunnel  No.  8,  using  roughened  calorimeters,  form  tlie  data 
base  for  the  rough  wall  transition  correlation  and  analysis.  In  Series  A, 
the  surface  roughness  varied  over  the  range  (0.6  to  80  mils)  so  that  the 
relative  roughness,  k/e,  spanned  the  values  of  interest  for  nosetip 
materials  in  flight.  Also,  values  of  the  wall  cooling  ratio  (T^/T^) 
and  momontum  thickness  Reynolds  number  covered  the  range  of  primary 
interest.  In  Series  J,  several  additional  calorimeter  models  of  various 
geometries  were  tested  in  order  to  evaluate  nose  radius  and  shape  effects 
on  transition;  all  models  had  the  same  surface  roughness  (3.5  mils). 
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Tha  data  froai  larlas  A wara  aorralatad  in  tarma  of  tha  noatantum 
thioKnaaa  Aaynolda  nuabar  (rapraaantlng  tha  ability  of  tha  boundary 
layar  to  amplify  diaturbaneaa)  varaua  a nodifiad  ralativa  roughnaaa, 
(Ta/T^)  (V*) « (rapraaantlng  tha  atrangth  of  tha  dlaturbanoa) . Thia 
oorralation  approach  la  llluatratad  on  tha  akateh  balowi 


two  aapaota  of  tranaition  wara  oorralatadi 

a Tranaition  onaat,  daaoribing  whan  roughnaaa  inducad  tranai- 
tion oooura  in  tha  aubaonio  ragion. 

a Tranaition  location > daaoribing  tha  poaition  of  tha  baginning 
on  tranaition  whan  it  oooura  in  tha  aubaonio  ragion. 

Thaaa  oritaria  may  ba  atatad  in  tarma  of  oritical  valuai  of  tha  tranai- 
tion paranatar, 
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M 215  , location  f 


(6-1) 


Sariaa  J data  aubatantlatad  both  criteria  davalopad  from  Sariaa  A data. 

An  analyaia  of  exact  amooth  wall  boundary  layer  profllaa  ahowad 
that  tha  ampirioal  dlaturbanoa  parameter,  (T./T„) (k/6) , correlataa  vary 
wall  with  tha  ralativa  kinetic  energy  at  tha  top  of  tha  roughnaaa  ala- 
manta,  along  with  phyaioal  raaaoning,  auggaata  that 

the  latter  paramatar  ia  tha  proper  phyaioal  maaauro  of  tha  dlaturbanoa. 
Uaing  thia  ooneapt,  a function. 


‘Evaluated  at  tha  aonic  point. 


(6-2) 
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was  davalopad  to  raplaca  in  tha  anpiirloal  disturbance  paramatar. 

thereby  accounting  for  tha  affects  of  blowing,  ablation,  wall  cooling 
ratio,  and/or  dissociation.  For  low  enthalpy,  nonablating  situations, 
tha  function  i|«  is  tha  wall  cooling  ratio  (T^T^) . For  the  high  enthalpy, 
ablating  case.  Equation  (6-2)  indicat<^s  the  affect  of  B*  and  P,/Py  on 
tha  kinetic  energy  profile  and,  hence,  on  the  roughnees  generated 
boundary  layer  disturbances.  The  transition  criteria  become 


( k V*’ 

( 255*, 

onset  ) 

Vipy)  " 

1 215  , 

location) 

Because  the  above  transition  oriteris  are  closely  associated 
with  physical  roughness  dimensions,  accurate  values  of  surface  roughness 
in  the  particular  flight  or  teat  being  analysed  are  essential.  In  clear 
air  flight  eases  and  some  ground  tests,  this  roughnuss  would  correspond 
to  that  associated  with  laminar  ablation.  In  tests  where  a nosetip  is 
rapidly  exposed  to  a severe  environment  or  is  exposed  to  particle  laden 
flows,  other  roughnesses  such  as  msohining  msrics,  fabrication  variations, 
or  impact  craters  may  be  important.  In  the  past,  clear  air  predictions 
have  been  based  on  affective  roughnesses  inferred  from  shape  change 
data.  Effective  roughness  dimensions  inferred  using  Equation  (6-3) 
should  agree  with  physical  roughness  dimensions  measured  from  photomicro- 
graphs of  laminar  ablation  models  if  roughness  does  not  change  during 
model  cooldown.  ^ 

For  ATJ-S  graphite,  the  rwghness  inferred  using  the  criteria 
given  by  Equation  (6-3)  above  is  between  0.4  and  0.6  mil  (Reference  6-1) 
ATJ-8  photomicrograph  data  are  reported  in  Reference  6-2  and  shown  in 
Figure  6-1.  Tha  two  lines  represent  the  best  fit  of  meaaurements  from 
high  (50  to  170  atm)  and  low  (0.05  to  25  atm)  stagnation  pressure  tost 
models.  High  pressure  tests  were  most  probably  turbulent;  whereas,  low 
pressure  tests  were  probably  laminar.  A pea)(-to-valley  roughness  height 
for  use  in  the  PANT  correlations  can  be  derived  from  Figure  6-1.  This 
roughness  is  dsfined  as  ths  dlstanco  from  the  significant  pealts  to  the 
significant  valleys  and  could  bo  represented  as  the  distance  between  the 
0.16  probability  datum  and  the  0.84  probability  da^um.  This  would, 
therefore,  include  the  middle  67  percent  of  all  peaks  or  valleys.  Ths 
peak-to-valley  roughness  (Kp^)  is 

kpv  ■ yj^g  - Ygg  ■ 0.67  - 0.18  - 0.5  mil  (laminar) 
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Figure  6-1.  Microroughneee  Height  Dietributlone  for  ATJ-S  Graphite 


For  an  upper  limit  (10  percent  to  90  percent) , 


k ■ 


^10  “ ^90  “ (upper) 


For  an  eatlmate  of  turbulent  microroughneee,  the  high  preeaure  curve 
glvea 

- 0.65  - 0.25  “ 0.6  mile  (turbulent) 


It  ie  encouraging  that  the  roughneae  height  inferred  from  the  ablation 
teat  data  (0.4  to  0.6  mil)  agreea  with  the  reaulte  of  the  photomicrograph 
meaaurementa . An  ATJ-S  roughneae  height  of  0.5  mil  la  therefore  recom- 
mended for  uae  with  the  newly  developed  correlation  in  the  abaence  of 
data  on  a apeciflc  material  aample. 

It  cannot  be  overemphaeized,  however,  that  for  each  flight  or 
ground  teat  of  graphite,  a aimilarly-prepared  apecimen,  which  haa  been 
aubjected  to  laminar  ablation  ahould  be  uaed  to  determine  the  apecific 
roughneaa  for  that  particular  caae.  An  accurate  roughneae  height 
meaaurement  la  particularly  important  becauae  noaetip  tranaltion  onset 
la  quite  aenaitive  to  laminar  roughness.  Figure  6-2  Illustrates  this 


Figure  6-2.  Treneition  Altitude  «s  a Function  of  Roughnaea  Height 
for  a Typical  XCBM  Trajectory 

aensltivlty  for  a typical  ICBM  trajectory  and  noaetlp  geometry.  A 
factor  of  2 decreaee  in  roughneaa  will  decrease  the  transition  altitude 
by  approximately  15 « 000  feet.  In  general « this  would  significantly 
reduce  nosetip  recession. 

6.2  IMPROVED  MODELING  OF  TURBULENT  HEAT  TRANSFER 

Transitional  and  turbulent  Stanton  numbers  are  coeiputed  using 
Reynolds  analogy  applied  to  the  results  of  streamwise  integrations  of 
the  integral  momentum  equation.  For  rough  surfaces,  these  heating 
rates  are  modified  to  account  for  roughness  augmentation.  Ablation 
effects  are  decoupled;  that  is,  blowing  corrections  are  computed  after 
the  nonablating  solution  huj  bean  generated. 

Two  aspects  of  the  technique  used  in  the  computation  of  turbu- 
lent heat  transfer  are  discussed  in  this  section.  In  Section  6.2.1, 
implementation  of  the  Crowell  composite  heating  model  (Reference  6-3) 
is  described.  This  model  provides  a realistic  low  Reynolds  number 
limit  to  the  turbulent  heat  transfer.  Section  6.2.2  treats  the  transi- 
tional heating  model.  In  Section  6.2.3,  rough  wall,  fully-turbulent 
heat  transfer  predictions  are  compared  to  rough  wall  calorimeter  data. 
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6.2.1  Crowll  Compesit*  H«*tinq  Modifiotion 

Previous  to  this  moAificatlonf  turbulsnt  hast  transfer  was  pra- 
dlcted  to  be  zero  at  the  stagnation  point  and  leas  than  laminar  heat 
transfer  for  momentum  thioknesa  Reynolds  numbers  below  about  50.  Thue, 
for  transition  at  low  Reynolds  numbers,  turbulent  heat  transfer  pre- 
dictions wore  physioally  unrealistic,  in  Reference  6-3,  Crowell  de- 
veloped a composite  skin  friction  expression  which,  when  used  in  the 
streamwise  integration  of  the  integral  stomentum  equation  (and  with 
Reynolds  analogy) , eliminated  the  low  Reynolds  number  anomaly.  The 
Reference  6-3  treatment  was  restricted  to  the  constant  property  smooth 
wall  case.  It  is  extended  herein  to  include  the  effects  of  variable 
properties  and  augmentation  of  heating  caused  by  surface  roughness. 

The  integral  momentum  equation  may  be  written  as 
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IPe«e®^ 
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(1  * H) 


(6-4) 


With  edge  conditions,  skin  friction  (t^) , and  boundary  layer  shape 
friction  (H)  defined  in  terms  of  the  other  variables,  Equation  (6-4) 
can  be  integrated  to  determine  the  momentum  thickness  as  a function  of 
streamwise  length,  6(s). 

Crowells  composite  expression  for  constant  property  smooth  wall 
skin  friction  is 


■'w  _ 0.222  . ( **e  \ 0.0128 


(6-5) 


This  is  shown  in  Figure  6-3  compared  to  the  laminar  and  turbulent  ex- 
pressions to  which  it  is  asymptotic.  For  RSg  ^ 200,  the  composite 
expression  is  within  1 percent  of  the  turbulent  skin  friction.  However, 
the  contribution  of  the  first  (laminar)  term  does  not  decrease  as 
rapidly  as  it  probably  should,  still  being  10  percent  at  Reg  - 1000. 

In  the  formulation  used  in  nosetip  analyses  turbulent  variable 
property  effects  are  modeled  via  Eckert  reference  enthalpy  properties 
(see  Reference  6-4) . The  reference  enthalpy  is  defined  as 


K - 0.5  h^  + 0.3  h,  + 0.2  hjj 


(6-6) 
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Figure  6>3.  Compoeite  Skin  Friction 
■o  that  the  reference  propertiee  are 

p(F,  p)  and  IT(K,  p) 


with 

pu.e 


In  the  dimensional  form  used  in  the  momentum  equation,  the  variable 
property  composite  skin  friction  is  written  as 


where 


0.222  U, 

— s — 


100  + Rig 
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Th«  r«f*r«no«  viMoaity  is  not  usod  in  th«  laminar  tarn  baeauaa  uaa  of 
tha  adga  viaooaity  yialda  a atoro  aseurata  momantum  thloknaaa  and  haat 
tranafar  pradlotion  at  tha  atagnation  jpoint. 

Tha  ahapa  factor,  H - 8*/e,  ramalna  to  ba  apaoifiad  bafora  solu- 
tion of  tha  intagral  momantum  aquatlona  ean  proeaad.  For  wall  cooling 
ratios  of  graataat  intaraat,  0.4  < T^T^  < 0.6,  tha  ahapa  factor  la  in 
tha  ranga,  0 < H < 1.0.  Tha  solution  of  tha  Intagral  momantum  aquation 
ia  not  particularly  sansitiva  to  variations  in  tha  ahapa  factor  and 
previous  versions  of  tha  ooda  used  H ■ in  order  to  allow  closed  form 
integration.  Xn  addition  to  its  direct  influence  In  the  integration  of 
tha  momantum  aquation,  H affaots  tha  value  of  tha  composite  momentum 
thioknasa  at  tha  stagnation  point,  tha  oomposita  momantum  aquation  at 
tha  stagnation  point  reduces  to 

6^(0)  - J ^(3  1 mp^Ju,7Jsy 

It  ia  not  posalbla  to  impose  an  initial  condition,  othar  than  this,  on 
tha  momantum  aquation.  Idaally,  0^(0)  should  ba  tha  same  as  tha  laminar 
moinantum  thickness  oaloulatad  by  the  ooda,  whloh  is 


" i Vp,Mu*/<!aT 

Tha  ratio  of  tha  two  stagnation  point  momantum  thicknesses  is 


(6-10) 


In  order  to  have  a raalistlc  initial  condition  for  tha  momantum  thlcknaas 
and  in  kaaping  with  previous  versions  of  the  coda,  a constant  ahapa 
factor,  H - -1,  is  assumed  in  tha  Integration  of  the  composite  momantum 
aquation. 
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Integration  of  the  momentum  equation  using  the  above  provides  a fully 
turbulent,  smooth  wall  solution  for  the  momentum  thickness  distribution  around 
the  body.  The  composite  smooth  wall  akin  friction  is  also  computed  during  the 
integration  process,  so  that  application  of  Reynolds  analogy  will  determine  the 
smooth  wall  heat  transfer. 

A laminar  Reynolds  analogy  factor,  Rj^,  is  applied  to  the  laminar  contri- 
bution to  the  composite  skin  friction  (Equation  (6-7)).  This  factor  is  used  to 
match  the  composite  heat  transfer  at  the  stagnation  point  to  the  Pay  and  Riddell 
(Reference  6-S)  prediction  used  in  the  laminar  part  of  the  code,  so  that 


— trriSaw'  — 


(6-11) 


The  turbulent  Reynolds  analogy  faotor,  R^,  is  ourrently  taken  to  be  unity 
over  the  entire  body.  There  is,  however,  some  evidence  from  turbulent  boundary 
layer  solutions*  that  R^  is  a function  of  pressure  gradient,  being  about  0.95 
on  the  nose  and  about  1.15  on  the  cone. 


The  smooth  wall  composite  heat  transfer  becomes 


PeVH,C,8  • \ 
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Surface  roughness,  if  sufficiently  large,  increases  laminar  as  well  as 
turbulent  heat  transfer.  These  effects  are  modeled  using  the  augmentation  fac- 
tors developed  in  Reference  6-6.  The  laminar  augmentation  factor,  F^,  is  the 
ratio  of  laminar  rough  surface  heat  transfer  to  laminar  smooth  surface  heat 
transfer,  both  for  the  same  environmental  and  surface  temperature  conditions. 
The  turbulent  augmentation  factor,  F^,  ie  the  ratio  of  heat  transfer  to  rough 
and  smooth  walls  in  fully  developed  turbulent  flow.  In  order  to  insure  that 
the  turbulent  augmentation  factor  becomes  dominant  as  the  composite  skin  fric- 
tion approaches  turbulent  skin  friction  the  augmentation  factors  are  applied 
as  shown  below. 


‘’e“e‘^H,C,R 


0.222u^ 
5 


100  \ 

100  . Si,  j * W.,t 


(6-13) 


^Boundary  Layer  integral  Matrix  Procedure  (BLIMP)  code. 
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If  th«  augmantatlon  faotora  wara  appliad  in  tha  mora  obvloua  way.  tha  Influanoa 
of  tha  laminar  augnantation  factor  would  paraiat  paat  tha  Ra^  whara  tha  flow 
would  bo  turbulant.  Thia  ia  aapaoially  aignifioant  for  vary  larga  roughnoaa 
and  for  waathar  oroaion  oaloulationa. 

For  vary  larga  roughnaaaaa  the  rough  wall  oompoaita  haat  tranafar  aa 
computed  from  Equation  (€-13)  falla  below  tha  rough  wall  laminar  prediction  in 
a amall  region  near  tha  atagnation  point.  Nhan  thia  oocurat  the  compoaita  haat 
tranafar  ia  aot  equal  to  tha  laminar  haat  tranafar. 


6.2.3  Tranaitional  Heating  Modal 

In  the  region  downatraom  of  tha  tranaition  points  tha  haat  transfer  is 
oaloulatod  aooording  to  tha  Parsh  transitional  heating  modal  (sea  Rafaranoa  6-7) 
Tha  transitional  Stanton  number  la  axprassad  as 


-H,T,R  ■ ^,C,R  “ 


(6-14) 


The  constant*  A*  ia  chosen  so  that  tha  Stanton  number  is  continuous. 

Thia  transitional  heating  modal  ia  aomawhat  simpler  than  tha  modal  which 
was  in  tha  coda  prior  to  tha  currant  effort.  As  will  ba  seen  from  tha  examples 
in  tha  next  section*  tha  present  modal  yields  quito  reasonable  pradiotlons  for 
rough  wall  ealorlmatars . These  predictions  are  batter  than  those  made  with  tha 
previous*  mora  complicated,  modal.  Further  rafinamants  would  ba  justifiable 
only  in  conjunction  with  an  improved  understanding  of  rough  surface  heating 
augmentation. 


6.2.3  Comparison  of  tha  Composite  and  Transitional  Predictions  to  Rough  Wall 
calorimeter  Bata  ^ * 

Rough  wall  calorimeter  data  from  Series  A and  Sarias  J taste  are  shown 
cosqparad  to  pradiotlons  in  Figures  6-4a  through  6-4f.  Thssa  figures  include 
examples  for  roughness  heights  of  0.6,  3.5*  and  40.0  milsi  nose  radii  of  1.5* 
2.5*  and  3.5  inohasi  with  wall  cooling  ratios  in  tha  range  0.43  < T^T^  < 0.74. 

Figures  6-4a  and  6-4b  show  t«ro  axaaiplas  from  Series  A with  small  and 
larga  roughnesses.  For  Run  105,  transition  is  pradictad  slightly  aft  of  its 
actual  position  which  results  in  a slight  undarpradiction  of  paa)c  heating. 

For  Run  174*  tha  actual  roughness  augmentation  of  paak  turbulant  heating  is 
slightly  larger  than  pradictad.  Soma  ralaminarisation  (net  medalad  in  tha 
present  coda)  is  evident  in  tha  region  of  tha  cone  tangent  point.  Both  Sorias 
examples  show  an  accurate  prediction  of  turbulant  heating  on  tha  cone. 
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Th«  n«)tt  four  tlguroi  show  sxsmplsa  frosi  Sories  J,  all  for  a roughnaaa 
height  of  3.S  mils.  Flguraa  6-4o  and  6-4d  ara  for  the  3.5-lnoh  nose  radius 
modal.  The  ovsrpradlotion  of  transitional  heating  for  Run  607  la  not  unique 
to  this  oalorlBstar.  This  was  obaarvsd  on  sevaral  runs  for  which  transition 
ooeurrsd  downstream  of  s/R^  > 0.4.  ividsntly,  the  same  forces  which  promote 
relaminiariaation  on  the  aft  spherical  surf sea  inhibit  turbulence  growth.  Run 
610  is  included  as  an  example  of  the  runs  for  which  the  total  temperature  waa 
low.  Agreement  of  the  prediction  with  the  data  is  quite  good.  Roth  Runs  607 
and  610  show  significant  relaminarisation. 

Pigurea  6-4e  and  6-4f  show  data  and  predictions  for  the  1.5-inch  nose 
radius  model.  These  examples  are  for  Run  620  at  test  times  of  1 and  10  aeconda. 
At  1 secondi  transition  is  prediotad  slightly  ahead  of  its  actual  positlcn 
loading  to  a mild  overprsdiction  of  peak  heating.  At  10  saoonds.  transition  is 
accurately  prediotad  but  the  composite  model  undarpredicts  peak  heating.  This 
is  probably  an  indication  that  the  dspendenoa  of  the  turbulent  roughness  augmen- 
tation factors  upon  wall  cooling  ratio  is  somewhat  in  error.  Because  of  the 
amall  nose  radius#  significant  relaminarisation  is  present  in  spite  of  the 
relatively  high  unit  Raynolds  number  of  this  example. 

In  general#  transition  location  and  transitional  heating  are  quite  wall 
predicted  for  the  rough  wall  calorimeters.  Heating  distributions  for  transition 
positions  between  s/Rj^  ■ 0.4  and  the  sonic  point  are  not  well  predicted.  Model- 
ing of  relaminarisation  is  lacking. 

6.3  SHAPE  CMANQB  MECHANICS 

A critically  important  aspect  of  nosetip  ablative  response  modeling  is 
the  technique  used  to  numerically  describe  the  instantaneous  shape  of  the  nose- 
tip  at  any  computational  time  step.  The  shape  definition  is  most  difficult  In 
the  stagnation  region  of  shapes  which  exhibit  a sharpened#  biconic  or  slender 
cone  stagnation  region.  During  the  PANT  program  numerous  tschniquas  for  model- 
ing the  shape  have  been  studied)  but  no  totally  satisfactory  technique  haa  been 
identified  to  date.  In  this  section#  the  various  techniques  are  described. 

The  problem  is  to  repreaent  the  body  shape  (local  surface  coordinates# 
slope#  and  curvature)  using  information  at  finite  nodal  locations.  The  complex- 
ity of  the  interacting  phenomena  and  associated  equations  means  that  many  of 
the  important  response  parameters  at  the  nodal  locations  must  be  computed 
explicitly.  That  is#  information  from  the  previous  time  step  must  be  used  to 
compute  current  time  results.  Por  example,  local  ablation  rates  at  a particular 
time  are  applied  in  a»st  shape  change  codes  over  a finite  computation  time  step. 
Since  criteria  for  explicit  extrapolation  of  instantaneoua  results  are  difficult 
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to  ••tablish,  nuMrioal  InstabllitlM  and  oaeillatloni  can  result.  Furthermora. 
ths  nunbar  of  surface  nodal  points  which  can  aeonoaiically  be  used  in  a shape 
change  calculation  is  liaitad.  The  nuaerical  techniques  which  circusivent  the 
difficulties  of  shape  description  are  here  termed  shape  change  numerics. 

Table  6-1  summsriiea  the  varioua  nusMrioal  techniques  incorporated  in 
the  nosetip  design  oemputar  oodaa  sodified  or  generated  during  the  PANT  program. 
Only  minor  changes  to  the  isodeling  of  aurfaoa  coordinate  evaluation  have 
occurred  during  the  development  of  the  PANT  codas.  Major  modifications  have 
been  associated  with  the  evaluation  of  the  local  surface  angle  at  the  nodal 
coordinates.  This  aspect  of  the  modeling  is  quite  Important  because  local 
pressures  and,  hence,  heating  rates  are  sensitive  to  the  local  surface  inclin- 
ation. The  smoothing  procedures  incorporated  in  the  lAAMT  and  PAOAM  cos^puter 
codes  provide  analytically  smooth  surface  angle  distributions  and  eliminate 
coordinate  location  anomalies  which  develop  from  calculational  inaccuracies. 

Zn  the  newly  developed  BPOI  computer  code,  nodal  location  adjustment  is  incon- 
venient because  eurface  points  are  tied  closely  to  the  transient  in-depth  con- 
duction grid.  Thus,  nodal  point  angles  are  determined  ao  that  numerical  shape 
instabilities  are  suppressed.  The  technique  is  described  in  Reference  6-12. 

Another  area  of  shape  change  numerics  which  has  been  modified  during  the 
PANT  program  is  the  definition  of  body  shape  in  the  near  stagnation  region. 
During  the  PANT  program  nosetip  recession  was  observed  to  produce  nearly  sharp 
bioonic  shapes.  Modeling  of  such  shapes  required  consideration  of  stagnation 
region  shape  details  which  were  on  a smaller  dimensional  seals  than  could  be 
economically  treated.  Special  nose  geometry  modeling  procedures  were  developed 
to  evaluate  the  instantaneous  nose  radius  not  from  the  nodal  point  locations, 
but  from  shape  stability  conditions.  Xn  the  BROS  code,  the  nose  rsdius  for 
sharp  shapas  is  evaluated  ao  that  the  stagnation  point  recession  rate  equals 
the  axial  recession  rate  of  the  first  surface  point  which  experiences  turbulent 
flow. 

6.4  BASSI.ZUB  PREDICTION  COMPARISONS 

Current  techniques  for  calculating  the  performance  of  a reentry  vehicle 
nosetip  require  aaaumptions  regarding  various  modeling  techniques.  The  assump- 
tions for  which  the  predictions  are  sensitive  are  itemiaed  belowt 

e The  shock  shape  evaluation  procedure  affects  the  solution  through 
the  influence  of  vortical  layer  swallowing  on  the  turbulent  boundary 
layer. 

e Local  surface  pressure  modeling  directly  influences  the  heat  tranafer 
calculation. 
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• Th«  •■■uai«d  skin  friction  law  used  in  tha  boundary  layar  momantum 
Intagral  aquation  af facts  tha  haat  transfar  calculation. 

a AasuMsd  transition  modallng  including  uatarlal  assumptions  such  as 
tha  affactiva  surface  roughness  is  important. 

a Tha  corralation  or  oorralations  used  to  modal  tha  affaota  of  surface 
micro-  and  maoro-roughnasa  on  heat  transfar  significantly  affact 
tha  rasponss  prediction. 

a Current  taohniquas  employ  relatively  arbitrary  spaoif ication  of 
maoro-roughnasB  height  and  growth  rata. 

a Assumptions  regarding  hataroganaous  surface  reactions  and  surface 
sublimation  are  important  because  tha  bhannoohamistry  models  datac- 
mina  tha  modeling  of  surface  particulate  removal  that  is  required  to 
match  ablation  data. 

a The  definition  of  the  nosatip  shape  profile  (i.a.,  point  location 
and  local  slope)  affects  all  subsaqusnt  calculations. 

Tha  above  list  represents  thosa  aspects  of  nosatip  modeling  which  are  las« 
than  oomplataly  understood.  Various  thaoratloal  and  semi-empirical  taohniquas 
have  provided  valuable  modeling  guidanoai  however.  Comparisons  to  available 
flight  data  provide  the  ultimate  validation  of  the  analysis  prooedure.  Thsre- 
fora«  seven  flight  cases  have  been  selected  by  SAM80  for  evaluation  of  the 
PANT  developed  computer  code  predictive  techniques.  These  flight  teat  oases 
provide  a wide  range  of  reentry  environments  against  which  the  merit  of  model- 
ing assumptions  may  be  assessed.  Section  6.4.1  describes  the  baseline  modeling 
and  corresponding  flight  predictions  obtained  at  the  start  of  the  effort  re- 
ported herein,  while  Section  6.4.2  presents  similar  results  obtained  using 
tha  current  modeling  updates. 


6.4.1  SAAHT  Code  Predictions 

Tha  SAANT  computer  code  predictions  were  performed  including  tha  shock 
shape,  pressure,  and  skin  friction  correlations  documented  in  Reference  6-10. 
Also,  the  body  curve-fit,  smoothing  option  described  in  this  reference  was 
used  for  the  baseline  predictions.  Tha  PANT  rough  well  nosatip  trnnsitlon 
model  based  on  displacement  thickness  was  used;  that  is. 


«R(f) 

I 


2000  location 
2300  onset 
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8ino«  th«  truialtlon  crltarlon  li  roughnaas  dapandant,  tha  corraet  mlcro- 

roughnaaa  muat  ba  apaclflad  in  ordar  to  aoeurataly  pradiot  tranaltlon  onaat 

and  location.  For  tha  graphitic  matariala  of  intaraat,  tha  appropriata  affac-v 

tiva  ffliororoughnaaa  for  uaa  with  tha  abova  tranaition  criterion  ia  0.4  mil  j 

(Rafaranca  6-1) . A diffarant  tranaition  crltarion  waa  applied  for  loeationa 

aft  of  tha  noaatip  ahouldar.  Thia  cona  tranaition  criterion  la  daacribad  aat 

I 

o.ntM  ! 

Ra^  “ 27Se  * at  tranaition 

J 

where  { 

■ edge  Mach  nunber  I 

Both  turbulent  and  laminar  heating  are  augmented  due  to  wall  roughnaaa 
by  tha  taehnigua  daacribad  in  Rafaranca  6-6.  Zn  laminar  flow,  the  material 
mlcroroughnaaa  in  utiliaad  with  tha  heating  augmentation  correlation.  In  tur- 
bulent flow  regimaa,  it  haa  bean  obnervad  that  aeallopa  develop  on  homogenaoua  . 

matariala.  Tha  affective  macroroughnaaa  (k)  for  aeallopa  haa  bean  aiodalad  in  I 

tha  banalina  pradictionn  aa 

k - 0.93P  "•*”  inch 

where 

P.n  • atagnation  praaaura  (pai) 

*2 

Macroroughnaaa  growth  ia  not  modeled  in  the  baaalina  predletiona.  However, 
downatraam  of  tha  noaatip  ahouldar.  the  macroroughnaaa  la  reduced  to  tha  micro- 
roughnaaa  value. 

The  thariRochemical  ablation  reaponae  for  the  graphitic  materlala  of 
intereat  wan  ahown  in  Reference  6-1  to  ba  accurately  modeled  by  aaauming  equi- 
librium diffualon  controlled  oxidation  and  aublimation.  The  Important  apeciea 
generated  by  heterogeneous  surface  reactions  include  CN.  C2N.  CjNj.  CO.  and 
CO2.  MO.  and  N2O.  although  the  chemical  species  allowed  in  the  modeling  of 
I carbon/alr  reactions  are  not  restricted  to  thasa.  Ground  test  evaluations  using 

tha  equilibrium  assumption  have  demonstrated  that  graphite  ablation  la  accurately 
modeled  assuming  no  mechanical  mass  removal.  Conaeguantly . the  baseline  flight 
predictions  do  not  incorporate  any  particulate  mass  ramoval. 

The  trajectory  and  geometry  parameters  for  the  seven  cases  are  presented 
in  normalized  form  in  Table  6-2.  Stagnation  point  recession  and  shape  profile 
histories  are  given  in  Figures  6-5a  through  6-Sg. 
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TABLE  6-2 

TRAJECTORY  AND  GEOHETRY  PAtAMETERS 


Case 

*^*REF 

%^^REF 

V'^REF 

«; 

(mc) 

a 

0.7S 

1.00 

1.00 

12.2 

159 

D 

1.50 

1.3S 

1.04 

8.2 

146 

0.75 

0.77 

0.90 

16.0 

123 

B 

1.25 

1.46 

1.04 

8.7 

138 

2.04 

0.81 

1.05 

17.0 

73 

a 

0.25 

1.23 

0.95 

11.0 

94 

B 

0.50 

0.27 

1.13 

Z6.0** 

110 

*Time  from  100  kft  to  LOS. 

**LOS  at  30  kft,  all  others  at  Impact. 


Figure  6-5.  Conti. 


Figure  6-5.  Continued 


Figure  6-S.  Continued 


Figure  6-5.  Concluded 
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€.4.2  1108  Cod<  Fr#diotioni 

Th«  8808  (Broalon  8h«p«-Ch*nga)  ooda  waa  davalopad  and  doeuMntad 
undar  8AN80/MBB8  dlraotion  partially  throuQh  a aupplaaiantal  agraaaiant  to 
anothar  oontraet  affort.  Coda  atodifioatlona  daaerlbad  in  daotiena  €.1,  6.2, 
and  6.3  wara  inoludad.  Tha  ooda  ehanpaa  dafinad  undar  tha  aupplaaantal  taaka 
allewad  aodaling  of  tha  folloainQi 

a ContinuuB  aroaion  aaaa  loaa 

a Haating  augnantation  in  aroaiva  anviromanta 

a Shock  layar  affaota  on  hydrcuMtaor  partiola  braakup  and  damiaa. 

a laprovad  shook  ahapa  avaluatior  taohniquaa 

a Transient  hast  conduction  via  DuPort-Frankal  schsaa 

Of  thasa  ohangas  only  tha  shock  shapa  nodaling  would  affact  tha  elaar  air, 
staady-stata  ablation  solutions  for  tha  savan  basalina  oasas.  Tha  upgraded 
shook  shapa  modal  is  dasoribad  in  Rafaranoa  6-13. 

Tha  modal ing  changas  dascrlbad  in  Saotions  6.1,  €.2,  and  6.3,  which 
aftaot  olaar  air  raoasaion,  nay  ba  sunmarisad  as  follows! 

a Transition  onsat  and  location  as  affaetad  by  blowing  and  dissocia- 
tion (6.1) 

a Turbulent  heat  transfer  formulation  (6.2) 
a Shapa  change  numerics  (6.3) 

Note  that  to  assess  tha  affect  of  the  transition  modal  change > tha  EROS  solu- 
tions ware  obtained  using  tha  same  nicroroughnass  for  transition  that  was  used 
in  the  8AAHT  solutions  (0.4  mil).  Thus,  transition  altitude  changas  are 
directly  rslatad  to  tha  criterion  formulation  change. 

Tha  only  additional  difference  between  tha  SAANT  basalina  solutions  and 
tha  baseline  solutions  generated  using  tha  EROS  cods  was  associated  with  the 
salaction  of  tha  turbulent  maeroroughnass  (scallop  affact  on  heat  transfer) . 

Tha  Sarias  H data  presented  in  Section  3 showed  that  the  use  of  actual  scallop 
dimensions  in  the  sandgrain  heating  correlation  significantly  overpradictad 
data.  Xn  tha  absence  of  definitive  data  for  the  material  of  interest  at  flight 
conditions,  a single  turbulent  region  rouglmass  was  selected.  The  value  for 
maeroroughnass  (k^)  was  taken  to  be 

k^  “ O.JOl  inch 


P 

K' 

p- 


6-27 


For  • ttafnatlon  praaaura  condition  of  100  ata,  tha  praaaura  dapandant  modal 
uaad  in  tha  SiUUlT  pradiotiona  givaa 

> 0.003  inch 

Tha  eorraaponding  raduction  in  tha  rou9>  wall  turbulant  heating  factor  ia  con- 
aiatant  with  tha  Sariaa  H raault. 

hacaaaion  and  ahapa  profila  hiatoriaa  from  tha  BROS  coda  pradictiona 
ara  givan  in  Figura  6-6a  through  6-6g.  Tha  ShAMT  baaalina  racaaaion  pradictiona 
ara  alao  ahown  in  tha  raapactiva  figuraa.  Tabla  6*3  proaanta  a compariaon  of 
pradictad  noaatip  tranaition  onaat  altitudaa  uaing  aavaral  tranaition  criteria. 
Concluaiona  from  tha  baaalina  pradictiona  nay  ba  aunmariaad  aa  followai 

a Tranaition  altitude  ia  predicted  alightly  higher  uaing  tha  modified 
tranaition  criterion. 

a Tranaition  altitudaa  ara  pradictad  lower  than  data  from  caaaa  A,  B 
(alightly) t and  Ci  higher  than  data  from  caae  D. 

Turbulent  ragima  racaaaion  rate  ia  predicted  to  be  aomewhat  lowers 
conaiatent  with  tha  modified  turbulant  heating  modal  and  acallop 
effects  modeling. 

a Turbulant  ragima  raeassion  rates  ara  in  batter  agreement  with  data 
using  EROS  modeling. 

a stagnation  point  racaaaion  histories  ara  more  monotonic  and  realistic 
because  of  improved  noae  region  numerics. 

a EROS  coda  shape  profilas  do  not  exhibit  drastic  inflaction  points 
and  are  generally  more  uniform  with  time  than  tha  SAANT  results. 
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TABLE  6-3 

SENSITIVITY  TO  TRANSITION  CORRE 


Sonic  Point 
Transition 
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8ECTI0H  7 


C0NCLU8ICN4 

Th«  oonolualons  from  th«  afforta  aummarliad  in  thla  voluna  of 
tha  Intarlm  Baport  for  tha  Paaaiva  Moaatip  Toohnology  (PANT)  program 
ara  itamiaad  balowt 

1.  Tha  APFDL  50  MN  are  jat  ablation  teats  of  varioua  alaa  and 
ahapa  noaatip  modala  providad  uaaful  ahapa  ohanga  and  bound- 
ary layer  tranaitlon  data.  Tha  ahapa  change  raaponaa  of  tha 
modala  waa  oonaiatant  with  ahapa  change  avanta  obaarvad  in 
numaroua  low-tamparatura  ablator  (LTA) , wind  tunnel  taata. 

2.  In  tha  AFFDL  SO  MW  aro  taata,  initial  modal  conf igurationa 
algniflcantly  affected  tha  early  time  ahapa  change  and  ra- 
ocaalon;  however,  once  fully-turbulant  profilaa  aatablaihad, 
raoaaalon  rataa  aqualiaad. 

3.  Heat  tranafar  maaauramanta  from  the  noaatip  replica  calorim- 
eter, which  waa  taatad  in  PANT  wind  tunnel  Sariaa  H,  pro- 
vided uaaful  data  on  tha  affecta  of  surface  scallops  on 
nosetlp  raaponaa.  Analysis  of  the  data  indicated  tha  follow- 
ing! 

a Heat  transfer  rates  to  the  scalloped  surface  calorimeter 
ara  signlf ioantly  higher  than  the  corresponding  smooth 
surface  predictions. 

a Tha  measured  heat  transfer  rates  generally  agree  with 
heating  rates  inferred  from  the  ablation  rates  of  the 
replicated  LTA  model. 

a The  roughness  agumentation  factor  (ratio  of  actual  to 
smooth  wall)  does  not  scale  as  Indicated  by  the  PANT 
sandgraln  roughness  effects  correlation,  as  anticipated. 

a Tha  augmentation  factor  at  specific  nosetip  locations  is 
nearly  independent  of  free  stream  Reynolds  number.  Loca- 
tions with  larger  roughness  indicate  greater  augmentation. 
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• Th«  Boallop  aurfao*  on  th«  Sorias  H oalorlnstsr  produoos 
loss  hasting  augmantstlon  than  a sandgraln  surfaea  having 
a sinilar  roughnass  haight.  Currant  aquivalant-sandgrain 
modals  do  not  axplain  tha  diffaraneas. 

4.  Tha  affaots  of  noaetlp  configuration  variations  and  trajectory 
dapandant  flow  conditions  on  nosatip  shaps  change  ware  anumar- 
atod  in  tha  Bar las  I Low-Teraparatura-Ablator  (LTA)  wind  tunnel 
tests.  These  constant  Reynolds  number  tests  indicated  that 
tha  Reynolds  numbar  condition  for  tha  irregular  shape  fonnatlon 
varies  significantly  with  stream  total  temperature,  initial 
nose  radius,  and  surface  roughnass. 

5.  Zn  tha  Series  I tests,  tha  formation  of  irregular  shapes  was 
accompanied  by  a significant  increase  in  the  vibrational 
accelerations  associated  with  unstable  or  quasl^unsteady 
flow  phenomena. 

6.  Zn  tha  Series  J wind  tunnel  tents,  heat  transfer  data  from 
roughened  calorlmaters . indicated  the  affects  of  body  size 
and  shape  on  boundary  layer  transition.  Analysis  showed 
that  the  transition  data  agreed  with  the  PANT  correlation 
developed  earlier  during  tha  analysis  of  Series  A.  2.5-inch 
nose  radius,  sphere/cona  calorimatar  data. 

7.  Detailed  study  of  boundary  layer  profile  calculations  indi- 
cated that  the  boundary  disturbance  which  initiates  transi- 
tion correlates  with  the  relative  flow  kinetic  energy  at  the 
top  of  a roughness  element.  Using  this  observation,  the 
transition  criterion  was  analytically  extended  to  flight 
condition  by  including  the  dependence  on  ablation  and  species 
dissociation.  Allowing  for  material  surface  roughness  un- 
certainty. the  upgraded  transition  model  agrees  with  avail- 
able flight  nosetlp  transition  data. 

8.  The  turbulent  heat  transfer  model  in  the  PANT  codes  was  up- 
graded to  include  a unified  skin  friction  formulation  which 
successfully  eliminated  low  Reynolds  nximber  turbulent  heat 
transfer  prediction  anomolies. 

9.  Seven  flight  test  cases  were  analyzed  to  determine  the  effects 
of  recent  shape  change  modeling  on  calculated  shape  profiles 
and  recession.  Zt  is  concluded  that  the  transition,  heat 
transfer,  and  code  numerics  modeling  changes  improved  the 
prediction  capability,  but  the  prediction  of  the  irregular 
shape  regime  remains  uncertain. 
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LIST  OF  SYMBOLS 
dlmtnslonlitt  bloMing  rate 
Stanton  nunbtr 
haat  tranifar  coafflclant 
roughnoss  height 
Mach  number  or  molecular  weight 
pressure 
note  radius 

normal  distance  from  axis 

PeUek/We 

OeUeB/Me 

streamwlse  distance  from  the  stagnation  point 
temperature 

streamwlse  velocity  component 
normal  velocity  component 
entry  spped 

normal  distance  from  the  surface 

dimensionless  velocity  gradient  or  ballistic  coefficient 
entry  angle 
displacement  thickness 
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LIST  OF  SYWQLS  (Concludid) 
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Ltvy-Lttt  noriMl  coordinate 
moaiintiia  thlckntsi 
dyntfilc  vUcosIty 
kliMmatIc  viscosity 
Livy-Laas  straa«Mlsa  coordinate 
mass  dansity 

affactiva  mall  cooling  ratio 


Subscripts 

a 

k 

0 

w 

Suoarscrlots 

* 


boundary  layar  adga 

top  of  roughnass  alamant 

modal  stagnation  condition  or  Incomprassibla  n 

wall 

sonic  condition 
fluctuation  quantity 

avaraga 
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A.1  INTRODUCTION 


Rough  Mall  transition  studies  were  undertaken  as  part  of  the  PANT  program 
In  order  to  obtain  a better  understanding  of  boundary  layer  transition  on  graphitic 
nosetips  during  reentry.  Previous  uncertainties  In  predicting  Nhen  boundary  layer 
transition  occurs  ahead  of  the  sonic  point  resulted  In  a factor  of  two  uncertainty 
In  predicted  nosetip  recession.  Nosetip  shape  change  (and  total  recession)  Is  very 
sensitive  to  the  onset  of  transition  In  the  subsonic  nose  region  but  It  also  depends 
on  the  location  of  transition  In  this  region  at  lower  altitudes. 

This  appendix  Is  a summary  of  the  experimental  and  theoretical  program  carried 
out  to  address  the  problem  of  predicting  rough  wall  transition  in  the  subsonic  region 
of  a graphitic  nosetip. 

In  the  remainder  of  this  section,  the  flight  transition  environment  and  Its 
simulation  In  ground  test  facilities  Is  discussed.  Section  A.2  summarizes  the  experi- 
mental program  which  was  comprised  of  two  hypersonic  wind  tunnel  test  series  utilizing 
rough,  nonporous,  thin  wall  calorimeters.  Transition  data  from  these  two  series  are 
listed  In  Tables  A-1  and  A-2.  Section  A. 3 treats  the  correlation  of  rough  wall  transi- 
tion data.  A modified  relative  roughness,  (T^/T^)(k/e) , Is  Identified  as  the  distur- 
bance parameter  which  correlates  with  the  momentum  thickness  Reynolds  number  at  transi- 
tion. A separate  correlation,  utilizing  the  same  parameters  evalauted  at  the  sonic 
point.  Is  presented  for  transition  onset,  that  Is,  the  condition  for  which  transition 
Is  Incipient  In  the  subsonic  region  of  a rough  nosetip.  Section  A. 4 Is  concerned  with 
a theoretical  extension  of  the  experimental  results  to  Include  the  effects  of  mass 
transfer  and  chemical  reactions.  The  empirically  established  disturbance  parameter, 

(T-/T„)(k/e),  Is  shown  to  correlate  with  the  relative  kinetic  energy  at  the  top  of 

o w 

the  roughness  elements,  P|(U|(VpgU^*,  which  Is  hypothesized  to  be  the  key  physical 
measure  of  the  disturbance.*  The  Influence  of  blowing  and  chemical  reaction  on  this 
physical  disturbance  parameter,  to  Infer  a modification  of 

the  empirical  disturbance  parameter,  (T^/T^)(k/e).  to  account  for  these  effects.  This 
modified  empirical  disturbance  parameter  replaces  (T^/T^)(k/e)  In  the  transition  cri- 
terion. 


A.  1.1  Flight  Transition  Environment 

Surface  roughness  may  be  expected  to  Influence  boundary  layer  transition  unless 
the  roughness  height  Is  small  In  comparison  to  the  boundary  layer  thickness.  Current 


*The  ratio  p|(Uk‘/PgU.’  Is  the  smooth  wall  kinetic  energy  at  a distance  of  one  roughness 
height  from  the  wall  divided  by  the  kinetic  energy  at  the  boundary  layer  edge.  The 
empirical  disturbance  parameter,  (T«/TM)(k/e),  Is  retained  In  the  prediction  procedure 
since  knowledge  of  the  values  p|(  and  u^  requires  detailed  boundary  layer  profiles  which 
are  not  available  In  the  shape  change  codes. 
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nosctlp  graphites  have  peak  to  valley  surface  roughnesses  in  the  range  0.1  mil  < k 
< 1.0  mil.  For  comparison  with  these  roughnesses,  the  laminar  sonic  point  momentum 
thickness  distribution  with  altitude  Is  shown  In  Figure  A-1  for  a 0.75*1nch  sphere- 
cone  nosetip  during  a typical  ballistic  reentry.  Since  the  roughness  Is  of  the  same 
order  of  magnitude  as  tho  momantian  thickness,  the  Influence  of  roughness  must  be 
taken  Into  account  when  predicting  boundary  Idyer  transition  on  graphite  nosetips. 

The  momentum  thickness  Reynolds  number,  Re^  ■ the  parameter  usually 

considered  to  reflect  the  tendency  of  disturbances  In  a boundary  layer  to  be  ampli- 
fied Into  turbulence.  This  parameter,  evaluated  at  the  sonic  point.  Is  shown  as  a 
function  of  altitude  In  Figure  A-1.  If  the  nosetip  had  a perfectly  smooth  wall, 
boundary  layer  transition  would  be  expected  to  occur  ahead  of  the  sonic  point  for 
ROg*  ~ 300.  Although  this  level  Is  not  attained  on  the  “typical”  0. 75-Inch  nose 
radius  reentry  body,  roughness  reduces  the  critical  value  uf  the  momentun  thickness 
Reynolds  number  required  for  transition.  The  results  detailed  In  later  sections 
Indicate  roughness  Induced  transition  on  this  nosetip  for  k > 0.1  mil. 

TWO  other  parameters  are  expected  to  Influence  nosetip  boundary  layer  transi- 
tion: the  wall  cooling  ratio,  T^/T^,  and  the  dimensionless  blowing  rate,  B'  ■ 

Sonic  point  values  of  these  parameters  are  shown  In  Figure  A-1.  The  Indicated  levels 
of  these  parameters  are  typical  of  graphite  ablation  In  flight. 

/'i.1.2  Previous  Data  Base 

Prior  to  the  PANT  program  most  sources  of  data  on  roughness  Induced  transition 
were  restricted  to  two-dimensional  adiabatic  flows  with  Isolated  roughness  elements. 
These  data  sources  are  described  by  Anderson  and  Bartlett  In  a survey  of  the  state 
of  the  art  at  the  beginning  of  the  PANT  program  (Reference  A-1).  There  were  a few 
sources  which  reported  calorimeter  data  for  the  situation  of  primary  Interest,  that 
Is,  uniform  distributed  roughness  on  a blunt  axisymroetric  body.  Data  from  these 
sources  were  Included  In  the  analysis  reported  by  Anderson  (Reference  A-2). 

j The  paucity  of  pertinent  data  Indicated  a need  for  a comprehensive  experimental 

) program  to  provide  an  adequate  collection  of  data  on  which  to  base  nosetip  transition 

criteria.  Ideally,  the  parameters  discussed  In  connection  with  Figure  A-1  (I.e., 
k/e,  Re^,  T^T^,  and  B'),  together  with  enthalpy  and  pressure  levels  typical  of 
reentry  would  be  closely  approximated.  Since  all  appropriate  parameter  variations 
cannot  be  simulated  simultaneously  In  any  ground  test  facility,  some  compromises  were 
necessary.  Three  types  of  facilities  were  considered  and  these  are  dlscuss'id  In  the 
following  section. 

A. 1.3  Candidate  Ground  Test  Facilities 

The  ballistic  range  offers  the  closest  simulation  of  the  flight  environment  as 
graphitic  models  can  be  flown  at  reentry  speeds  and  pressures.  Because  of  the  high 
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vtlocItiM  and  short  tost  tlms.  InitruMntatlon  Is  difficult;  hoMOvor,  optical 
ttchniquas  art  bolng  davolopad.  Modal  prohaatlng  and  modal  racovary  tachniquas 
ara  also  baing  davalopad. 

High  powar  (50  MU)  arcjats  provida  simulation  of  flight  prassura  levels  and 
fairly  high  enthalpy  levels  In  which  to  test  graphite  rnodals.  Tests  may  ba  run  for 
several  seconds  at  constant  conditions.  The  primary  shortcomings  of  this  type  of 
facility  are  the  nonunifonn  enthalpy  and  entropy  distributions  In  the  stream  and 
the  high  level  of  free>stream  turbulence  and  contamination. 

The  hypersonic  wind  tunnel,  although  It  operates  at  low  pressure  and  enthalpy 
In  comparison  to  flight,  allows  simulation  of  several  of  the  most  Important  dimen- 
sionless transition  parameters.  Using  nonporous,  thin  wall,  metal  calorimeters  with 
roughness  heights  of  the  order  of  1 to  10  mils,  the  ranges  of  Interest  of  the  momen- 
tum thickness  Reynolds  number,  the  ratio  of  roughness  height  to  momentum  thickness 
(relative  roughness)  and  the  wall  cooling  ratio  can  all  be  simulated.  The  test 
stream  Is  clean  and  wall  defined.  Also,  Instrumentation  Is  relatively  straightforward 
and  test  times  are  conveniently  long.  For  tiMsa  reasons  the  hypersonic  wind  tunnel 
was  chosen  for  the  rough  wall  transition  experiments.  The  effects  of  high  enthalpy 
level  end  ablation  are  not  simulated;  however,  as  discussed  In  Section  A.4,  these 
effects  may  be  treated  theoretically. 

After  a survey  of  available  wind  tunnels,  NOL  tunnel  number  8 was  selected 
as  the  facility  which  provided  the  best  transition  testing  environment  (see  Refer- 
ence A-3).  Operating  at  Mach  5,  the  tunnel  unit  Reynolds  number  range  was  more  than 
adequate  to  ensure  the  capability  of  obtaining  both  laminar  and  turbulent  flows  In 
the  subsonic  region  of  models  with  roughnesses  In  the  range  of  0.6  to  40  mils.  At 
the  same  time,  the  momentum  thicknesses  were  such  that  the  relative  roughness  , k/e, 
covered  the  range  of  Interest  for  flight.  The  tunnel  may  be  run  at  stagnation  tem- 
perature up  to  at  least  1400*R  so  that  with  the  model  Initially  at  room  temperature, 
wall  cooling  ratios  In  the  range  of  0.4  to  0.8  are  experienced  as  the  model  Is 
heated.  In  summary,  this  facility  allows  the  close  approximation  of  flight  values 
of  the  dimensionless  parameters  which  have  the  greatest  Influence  on  transition. 

A. 2 EXPERIMENTAL  PROGRAM 

Two  series  of  rough  wall  calorimeter  transition  experiments  ware  executed. 

The  first  r:  these.  Series  A,  had  a second  objective,  namely  to  determine  heating 
augmentation  due  to  roughness  (see  Reference  A-4).  In  that  series  several  sphere 
cones  and  one  blunt  cone  were  tested.  All  were  nominally  the  same  size  with  a wide 
variation  In  roughness  height.  In  the  second  series.  Series  J,  several  new  models 
of  various  sizes  and  shapes  were  utilized.  All  of  these  had  the  same  roughness  height. 


A.2.1  T«it  Str1«i  A 

A. 2. 1.1  Mbdal  Fabrication  and  Instrumantatlon 


Tht  modal  gaomitrlat  for  this  test  sarlas  art  shown  In  Figure  A>2.  Seven 
rough  wall  sphere  cones  were  tested,  all  with  ■ 2,5  Inches  and  cone  half  angles 

of  8*.  The  eighth  rough  wall  model  was  a blunt  cone,  nominally  the  same  size,  with 
■ 5.0  Inches,  a 0.5*1nch  shoulder  radius  and  a cone  half  angle  of  S'*.  The  models 
were  machined  from  elactroformed  nickel. 

Surface  roughnHses  on  four  models  were  generated  using  a grit  blasting  tech- 
nique; three  sphere  cones  with  peak  to  valley  roughnesses  of  0.6,  1.6,  and  3.0  mils 
and  the  blunt  cone  with  a roughness  of  3.0  mils.  The  remaining  four  sphere  cones 
had  copper  particles  furnace  brazed  to  their  surface  with  resulting  roughness  heights 
of  3.0,  10.0,  40.0  and  60.0  mils.  The  two  roughest  models  were  primarily  Intended 
to  provide  data  on  the  effects  of  roughness  on  turbulent  heat  transfer. 

Figure  A-3  shows  photomicrographs  of  a sample  3-m11  grit  blasted  surface  i 

and  a sample  3-m11  brazed  particle  surface.  The  peak  to  valley  roughness  height  | 

Is  the  distance  between  the  significant  peaks  and  the  significant  valleys,  but  Is 
not  the  distance  from  the  very  highest  peak  to  the  lowest  valley.  The  character  j 

of  an  ablated  graphite  surface  Is  probably  best  represented  by  the  grit  blasted 
surface. 

In  order  to  provide  an  evaluation  of  the  subjective  determination  of  roughness 
height  described  above,  samples  of  all  surfaces  with  k < 10  mil  were  analyzed  using 
a Talysurf  profllometer.*  The  RNS  roughnMs  determined  In  this  way  Is  shown  compared 
to  the  peak  to  valley  roughness  heights  In  Figure  A-4.  One  sample  from  Series  J Is 
Included.  Apparently  the  peak  to  valley  measure  Is  equal  to  four  times  the  RMS 
roughness. 

The  models  were  Instrumented  with  chromel-alumel  thermocouples  spot  welded  to 
the  Interior  surface.  TWo  or  three  primary  rays  on  each  model  had  thermocouples 
located  at  every  O.S-Inch  of  the  streamwise  dimension.  Other  thermocouples,  less 
closely  spaced,  were  placed  on  other  rays  In  order  to  assess  asymmetries  in  the  heat 
transfer  to  the  models. 

i Further  details  of  model  fabrication  and  Instrumentation  can  be  found  In 

Reference  A-2  or  Reference  A- 5. 

(• 

I 

A. 2. 1.2  Data  Reduction 

I 

i During  the  tests  thermocouple  outputs  were  recorded  at  the  rate  of  about  18 

I times  per  second.  These  temperature  data  were  temporally  smoothed  using  a quadratic 

*The  results  of  these  measurements  were  provided  by  Mr.  Gary  Denman  of  AFML/MXS. 
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splint  fit  proctdurt.  Afttr  an  Initial  rttponsa  tint  which  Is  Itss  than  0.5  sacondt 
front  and  back  wall  tanparatura  risa  ratas  art  assantlally  aqual.  Tharafora.  haat 
transfer  to  tha  calorlmatar  was  calculated  consldarlng  only  anargy  storage  associated 
with  the  tanparatura  risa  rata  and  straaiMlsa  conduction.  Details  of  tha  data  reduc- 
tion procedures  are  given  In  Rafamnca  A-2. 

For  tha  purpose  of  datarartnlng  transition  locations,  distributions  of  axparl- 
mental  haat  transfer  coefficients  were  plotted  at  several  test  times  for  comparison 
with  tha  laminar  distribution.  The  laminar  distribution  Is  nearly  Independent  of 
wall  temperature  and,  therefore,  test  time.  The  transition  location  was  defined  as 
tha  point  where  tha  experimental  distribution  departed  from  the  laminar  dlstrlbutlom 
see  Section  A.2.2.2  and  Figure  A-6.  Several  examplK  are  shown  In  Reference  A-2. 

Information  describing  the  test  conditions  for  the  runs  yielding  transition 
data  Is  presented  In  Table  A-1.  Transition  locations  at  three  (or  sometimes  two) 
test  times  along  with  other  pertinent  data  are  also  shown.  Correlations  of  exact  smooth 
wall  laminar  boundary  layer  solutions  were  used  to  calculate  the  boundary  layer  param- 
eters e and  d*.  The  details  of  these  methods  are  given  In  Reference  A-2. 

The  correlation  of  this  data  la  described  In  Section  A. 3.  In  order  to  deter- 
mine whether  any  Important  physics  associated  with  body  size  and  shape  were  lacking 
from  the  correlation,  a second  set  of  experiments.  Series  J,  described  below,  was 
performed. 

A. 2. 2 Test  Series  J 

A. 2. 2.1  Model  Fabrication  and  Instrumentation 

The  model  geometries  for  this  test  series  are  shown  In  Figure  A-5.  Three  new 
sphere  cone  models  were  fabricated.  These  had  nose  radii  of  0.75,  1.5,  and  3.5  Inches, 
all  with  an  a*  cone  half  angle.  The  two  other  models  were  a 60*  bIconIc  and  a laminar 
ablated  shape.  All  of  the  Series  J models  were  machined  from  cold  rolled  nickel  and 
ware  grit  blasted  to  a peak  to  valley  roughness  height  of  3.5  mils  (sea  Figure  A-4). 

Backfaca  chromel-alumel  thermocouples  were  again  used.  Thermocouple  locations 
varied  from  model  to  model  but  were  roughly  similar  to  the  Installation  on  the  Series  A 
models. 

For  more  detailed  Information  regarding  the  Series  0 models  and  Instrumentation, 
see  Reference  A-6. 

A.2.2.2  Data  Reduction 

Data  recording  and  smoothing  techniques  as  well  as  the  computation  of  heat 
transfer  coefficients  were  essentially  the  same  as  they  were  for  Series  A.  An  example 
plot  of  heat  transfer  coefficient  distributions  Is  shown  In  Figure  A-6.  Transition 


TABLE  A-1 

SERIES  A TRANSITION  LOCATION  DATA 
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2.60  1240  4.21 


2.01  1245  3.24 


2.5  10.0  0.39  1200 


Tim#  «/RnV'^#  U./V. 

(MC)  OOVft) 
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locations  for  1.  5.  and  10  second  test  timas  are  Indicated  by  the  Intersection  of  the 
straight  lines  with  the  laminar  distribution.  Run  conditions  and  transition  location 
data  for  Series  J are  listed  In  Table  A-2.  Data  from  the  0. 75-Inch  nose  radius  model 
were  not  used  here  but  were  used  In  the  onset  correlation  (see  Section  A. 3. 3), 

A. 3 ROUOKNESS  INDUCED  TRANSITION  CORRELATION 
A. 3.1  Modeling  of  Roughness  Induced  Transititw 

In  the  present  analysis  It  Is  assumed  that  the  response  of  a laminar  boundary 
layer  to  a disturbance  will  depend  primarily  upon  the  magnitude  of  the  disturbance 
and  the  potential  of  the  boundary  layer  to  amplify  disturbances.  With  proper  Identi- 
fication of  the  disturbance  and  the  amplification  potential,  the  boundary  layer  would 
be  expected  to  respond  as  Indicated  In  the  sketch.  Small  disturbances  would  require 


f 


! 

( 


t 


larger  amplification  potential  to  Initiate  transition  than  would  large  (ilsturbances. 

This  sort  of  Idea  was  used  to  correlate  transition  data  from  low  speed 
boundary  layers  on  adiabatic  flat  plates  which  are  subject  to  distrubance  from  the 
wind  tunnel  free  stream  turbulence.  The  data,  which  were  taken  from  Schllchting 
(Reference  A-7)  are  from  experiments  carried  out  by  several  Investigators.  These 
data  are  shown  In  Figure  A-7.  For  this  case,  the  free  stream  turbulence  Intensity 
Is  taken  to  be  the  measure  of  the  disturbance,  and  It  apparently  determines  the  value 
of  ROg  (amplification  potential)  at  which  transition  occurs. 

The  turbulence  Intensity  (abscissa  In  Figure  A-7)  Is  the  square  root  of  the 
relative  kinetic  energy  of  the  turbulent  fluctuations.  This  Is  perhaps  an  Indication 
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TABLE  A-2 

SERIES  J TRANSITIOH  LOCATION  DATA 
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1 
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MM 
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10 
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0.69 
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m 
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0.68 

1 
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5 
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10 
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1 
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10 
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1 
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0.50 

10 
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1.79 
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0.78 
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1 
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0.61 

1.05 
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10 
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m 
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0.43 
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(lOVft) 

1.5 

3.5  0.78 

975 

1.68 
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that.  In  ganaral,  the  klnatlc  tnargy  of  a dlsturbanca  Is  the  proper  Indicator  of  Its 
magnitude.  This  concept  will  be  used  In  a later  section  where  the  correlation  de- 
veloped below  Is  extended  theoretically  to  Include  blowing. 

A.3.1.1  Amplification  Potential  Representation 

Following  this  example,  the  momentun  thickness  Reynolds  number  (evaluated  as 
If  the  wall  were  smooth)  Is  used  as  the  parameter  whose  magnitude  represents  the 
tendency  of  the  boundary  layer  to  amplify  disturbances.  FOr  the  above  constant  pro- 
perty example,  the  displacement  thickness  could  have  been  used  Instead  of  the  momen- 
tum thickness  because  the  shape  factor  Is  constant  (6Ve  ■ 2.59).  However,  In  the 
subsonic  region  of  a blunt  body  the  displacement  thickness  Is  a strong  function  of 
wall  cooling  ratio  while  the  momentum  thickness  Is  relatively  Insensitive  to  that 
parameter  (see  Appendix  C,  Reference  A-2).  Anderson  (Reference  A-2)  discussed  corre- 
lations based  on  6*  along  with  those  based  on  e.  The  correlations  based  on  displace- 
ment thickness  are  no  longer  reconmepded  because  (a)  6*  < 0 for  T^T^  < 0.2,  (b)  e Is 
most  readily  obtainable  from  engineering  codes,  and  (c)  the  correlation  based  on  a 
momentum  thickness  Is  slightly  better  when  the  data  from  Series  J are  Included. 

In  the  following  analysis,  the  momentum  thickness  Reynolds  number  Is  used  as  the 
sole  measure  of  the  amplification  potential.  However,  two  other  Influences,  known 
to  be  stabilizing,  are  at  work.  These  are  the  favorable  streamwlse  pressure  gradient 
and  the  normal  pressure  gradient  associated  with  the  centrifugal  acceleration.  It 
was  not  found  necessary  to  explicitly  Include  either  of  these  effects  In  the  corre- 
lation. However,  the  difference  between  the  onset  and  location  criteria  discussed  In 
a later  section  1$  most  likely  the  result  of  these  effects. 

A. 3. 1.2  Disturbance  Model 

As  mentioned  above,  the  flat  plate  example  (Figure  A-7)  suggests  that  the 
kinetic  energy  of  a disturbance  Is  the  proper  measure  of  Its  strength  to  promote 
transition.  In  the  absence  of  a detailed  knowledge  of  the  flow  around  roughness 
elements  Inmersed  In  a shear  flow,  this  kinetic  energy  Is  Impossible  to  evaluate. 
Instead,  a correlation  Is  developed  In  terms  of  the  simple  parameters  available  from 
an  engineering  code  based  on  the  Integral  momentum  equation.  The  simple  disturbance 
parameter  Identified  during  this  correlation  will  later  be  shown  to  be  associated 
with  the  kinetic  energy  of  the  mean  boundary  layer  flow  at  the  top  of  the  roughness 
element. 

The  basic  disturbance  parameter  Is  taken  to  be  the  relative  roughness  defined 
as  k/e.  This  Is  physically  meaningful  because  roughnesses  which  protrude  further 
out  In  a boundary  layer  will  cause  relatively  larger  disturbances.  The  magnitude 
of  this  Influence  depends  upon  the  density  and  velocity  profiles  so  that  the 
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boundary  condltloni  and  8'  can  also  ba  axpaetad  to  havo  an  effect  on  the 
strength  of  the  disturbance.  The  role  of  T^/T^  Is  evident  from  the  data  while  that 
of  B*  Is  Inferred  theoretically. 

A. 3. 1.3  Transition  Parameter  Distributions  Around  a Sphere  Cone 

Distributions  of  momentum  thckness*  boundary  layer  edge  unit  Reynolds  number, 
well  cooling  ratio,  and  momentum  thIcknMS  Reynolds  number  around  a sphere  cone  at 
uniform  wall  temperature  are  shown  In  Figure  A-8.  ROg  Is  seen  to  be  a smoothly 
Increasing  function  of  streamwise  distance  because  of  compensating  trends  In  the 
unit  Reynolds  number  and  momentum  thickness  distributions. 

For  uniform  roughness  height,  the  growth  of  the  momentum  thickness  tends  to 
decrease  the  disturbance  (k/e)  while  Increasing  the  amplification  potential  (Reg). 

If  these  effects  were  exactly  compensating  (they  aren't  quite)  the  location  of  transi* 
tion  would  not  depend  on  the  momentum  thickness  at  all.  In  that  case.  Ignoring  wall 
cooling  ratio  effects,  the  correlation  would  be  of  the  form 
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that  Is,  a critical  roughness  Reynolds  number.  This  may  ba  regarded  as  a reasonable 
approximate  Indicator  of  the  variation  of  the  severity  of  the  transition  environment 
around  the  body*  The  distribution  of  Rej^,  which  Is  the  same  as  the  unit  Reynolds 
number.  Indicates  an  Increasing  severity  of  transition  environment  with  streaiiMlse 
distance  until  a maximum  near  the  sonic  point  followed  by  a 30  percent  reduction  to  the 
cone  level.  The  critical  Re|^  transition  model  would  therefore  Indicate  that  transi- 
tion should  only  occur  ahead  of  the  sonic  point.  PANT  data.  In  fact,  show  that  on  the 
spherical  surface  rough  wall  transition  locations  are  restricted  to  the  subsonic 
region. 

A.3.2  Transition  Location  Data  Correlation 

Initial  efforts  to  correlate  transition  location  data  were  In  terms  of  Re„ 
as  the  amplification  parameter  versus  k/e  as  the  disturbance  parameter.  This  Is 
shown  qualitatively  In  the  sketch  In  the  top  left  of  page  A-24.  The  primary  short- 
coming of  this  method  of  correlation  Is  that  the  data  trend  during  a run  as  the 
wall  temperature  Increases  Is  not  parallel  to  the  main  (run  to  run)  data  trend.  The 
simple  R6|^  criterion  mentioned  above  Is  also  shown  on  this  sketch.  It  does  not  have 
the  correct  slope  to  properly  correlate  the  data  at  any  particular  wall  cooling  ratio. 
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It  was  found  that  tha  wall  cooling  ratio  tlnaa  tha  momantum  thicknati  It  an 
affactiva  boundary  layar  thicknatt  which  allowad  a tatltfactory  data  corralatlon. 

This  schama  Is  shown  In  tha  skatch  on  tha  right.  Data  trands  for  various  runt  (nota 
larga  roughnMt  axcaptlon  balow)  ara  atsantlally  collnaar.  This  corralatlon  tachniqua 
Is  utad  throughout  what  follows. 

A.3.2.1  Corralatlon  of  Sarlat  A Transition  Location  Data 

Tha  Sarlas  A data  llstad  In  Tabla  A-1  Is  shown  on  tha  abova  suggastad  corra* 
latlon  plana  In  Figura  A-9.  In  Rafaranca  A>2»  data  from  two  rays  wars  usually  ra- 
portad;  hara.  fbr  tha  saka  of  bravlty*  only  ona  ray  It  utad.  Modal  gsomatrlas  for 
this  sarlas  ara  shown  In  Figura  A-2.  Thata  modal s had  roughnats  halghts  batwaan 
0.6  and  40.0  mils  as  Indicatad  In  Figura  A-9.  A11  of  tha  data  axcapt  that  from  tha 
40.0>m11  modal  ara  corralated  within  20  parcant  by  tha  aquation 


Tha  axprasslon  on  tha  1aft  of  this  aquation  will  haralnaftar  ba  callad  tha  transition 
parametar. 

Thara  It  soma  possibility  that  another  source  of  disturbance,  namely,  wind 
tunnel  noise.  1s  slightly  Influencing  tha  data  from  the  modal  with  tha  lowest  rough- 
ness. Transition  Reynolds  numbers  for  these  data  appear  to  ba  a little  low  with 
respect  to  the  correlation  line.  If  this  Is  Indeed  a noise  effect  It  would  Indicate 
that  tha  roughness  on  this  modal  Is  the  threshold  for  roughness  domination  for  tha 
noise  levels  ancounterad  In  this  tunnel. 


It  Is  not  too  surprising  that  tha  data  from  tha  40  mil  modal  do  not  agraa 
with  tha  corralatlon  sines  k/e  Is  In  tha  ranga  of  20  to  30,  so  that  tha  roughnass  Is 
largar  than  tha  smooth  wall  boundary  laytr  thicknass.  It  can  ba  saan  from  Fig* 
ura  A-1  that  thasa  ralativa  roughnassas  ars  much  largar  than  thosa  ancountarad  In 
flight  with  k < 1.0  mil  and,  tharafora,  thasa  data  ara  not  of  too  much  concarn. 

As  mantlonad  In  Sactlon  A.2.1,  It  was  Judgad  nacassary,  aftar  tha  corralatlon 
of  tha  Sarlas  A data,  to  enlarga  tha  data  basa  to  Includa  exparlmants  with  a widar 
ranga  of  nosatlp  slia  and  shapa.  For  this  purposa,  Tast  Sarlas  J was  parformad. 


A.3.2.1  Corralatlon  of  Sarlas  J Transition  Location  Dita 

Tha  corralatlon  of  tha  Sarlas  J transition  location  data  (llstad  In  Tabla  A-2) 
Is  shown  In  Flgura  A-10.  Modal  gaomatrlas  ara  shown  In  Figura  A*3.  Excapt  for  tha 
ona  Sarlas  A modal  ■ 2.S  Inch,  k ■ 3.0  mil)  which  was  rausad,  all  of  thasa  modals 
had  roughnassas  of  3.5  mil.  No  axplldt  Influanca  of  body  slza  and  shapa  Is  avidant, 
and  tha  quality  of  tha  corralatlon  Is  about  tha  same  at  It  It  for  tha  Sarlas  A data. 

Although  the  data  from  tha  bIconIc  modal  fall  within  the  scatter  band,  tha 
trend  with  Increasing  wall  tamparatura  Is  not  parallel  to  tha  correlation.  This  It 
not  a typical  nosatlp  shape,  at  least  not  until  aftar  transition  and  sharpening  have 
occurred.  Thasa  data  would  ba  useful  at  part  of  a cone  transition  data  bate. 


A. 3. 3 Correlation  of  Transition  Onset  Data 


Tha  correlation  presented  In  Sactlon  A.2  Indlcatas  tha  location  of  transition 
whan  It  occurs.  Howavar,  roughnass  dominated  transition  Is  not  observed  to  occur 
In  a region  Just  forward  of  the  sonic  point  even  when  the  transition  parameter 
exceeds  215.  This  anomaly  Is  rationalized  through  the  use  of  a separata  onset  cri- 
terion describing  the  conditions  at  which  the  transition  location  Just  enters  (or 
loaves)  tha  subsonic  region.  As  mentioned  In  Sactlon  A.1,  transition  onset  Is  tha 
most  Important  single  event  In  determining  total  nosatlp  recession  and  therefore 
deserves  specific  attention. 

Flgura  A-11  shows  distributions  of  tha  transition  parameter  around  a sphere 
cone  for  two  conditions  which  differ  slightly  In  tha  severity  of  tha  transition 
anvtroiment.  These  are  computed  assuming  uniform  wall  temperature  and  uniform  rough* 
n«s*  lie'ght  Note  that  thosa  distributions  ara  quite  similar  to  distributions  of 
Ra|j.  For  tha  conditions  represented  by  tha  lower  curve,  transition  Is  not  observed 
(except  perhaps  on  the  cone)  even  though  the  transition  parameter  exceeds  215.  As 
tha  environment  becomes  more  severe,  say  by  Increasing  the  sonic  point  unit  Reynolds 
number,  transition  first  appears  whan  tha  transition  paramatar  evaluated  at  tha  sonic 
point  becomes  255.  Tha  location  criterion  than  determines  the  transition  point. 


Figure  A-10.  Correlation  of  Series  J transition  location  data 
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Distributions  of  the  transition  paraaeter  around  a sphara 
of  unifom  uall  taaperature. 


For  iny  singit  goomttry,  thi  onsot  critorlon  could  bo  replaced  by  a geometri- 
cal roftrictlon  on  the  application  of  the  location  criterion.  HoMver,  the  present 
scheme  appears  to  be  applicable  to  blunter  shapes  as  well;  for  example,  the  Series  A 
blunt  cone. 

Transition  onset  data  from  both  test  series  are  shown  In  Figures  A-12.  Data 
from  the  bIconIc  model  and  the  sphere  cone  with  a 40  mil  roughness  height  are  not 
Included.  The  various  runs  are  classified  according  to  whether  transition  was  located 
In  the  subsonic  region  during  the  entire  run  (turbulent),  transition  passes  aft  of  the 
sonic  point  during  a run  (onset),  or,  no  subsonic  transition  was  observed  (laminar). 
For  laminar  runs,  sonic  values  of  the  parameters  are  evaluated  at  the  beginning  of  the 
run  when  conditions  are  moat  severe.  For  onset  runs,  the  parameters  are  evalauted  at 
onset  time  and  for  turbulent  runs  they  ere  evalauted  at  the  latest  data  time;  that  Is, 
when  conditions  are  least  severe  (largest  T^T^). 

In  general,  the  onset  data  are  extremely  well  correlated  by  the  suggested 
expression.  Minor  exceptions  are  apparent  for  the  model  with  the  smallest  roughness. 
Again,  since  the  Reynolds  numbers  for  those  data  are  low,  this  may  be  an  Indication 
of  a minor  tunnel  noise  Influence  at  this  small  relative  roughness. 


A.  3. 4 Sunaaryoftorrel  atl  on  Res  u1  ts 


Using  Re^  as  the  amplification  parameter,  a modified  relative  roughness, 
(T^/T^)(k/e),  has  been  empirically  Identified  as  the  disturbance  parameter.  The 
correlations  for  onset  and  location  may  be  stated  In  terms  of  critical  values  of  the 
transition  parameter 
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location 


The  data  correlated  were  from  the  subsonic  region  of  roughened  blunt  calorimeters 
tested  In  a hypersonic  wind  tunnel.  Tests  covered  the  ranges  of  pertinent  parameters 
of  Interest  for  graphite  In  flight  with  the  exception  of  pressure,  enthalpy  level, 
and  blowing  rate. 

A.4  THEORETICAL  EXTENSION  OF  DISTURBANCE  MODELING  TO  FLIGHT  CONDITIONS 

In  the  previous  section,  the  modified  relative  roughness.  (T./T„)(k/e) . was 
empirically  established  as  the  measure  of  the  disturbance  Introduced  by  roughness 
elements  on  solid  calorimeters.  This  disturbance  parameter  Is  believed  to  be  repre- 
sentative of  a more  physically  meaningful  measure  of  the  disturbance,  such  as  the 


^Evaluated  at  the  sonic  point. 


Transition  onset  data  correlation 


rtlativ*  kinttic  entrgy  of  thi  disturbance  (refer  to  the  flat  plate  example  of  Figure 
A-7).  By  analogy  with  this  example,  and  with  the  additional  postulate  that  the  kinetic 
energy  of  the  disturbance  Introduced  by  the  roughness  Is  closely  related  to  the  mean 
kinetic  energy  at  the  top  of  the  roughness  elements,  the  following  |‘(ypothes1t  Is 

suggested: 

The  relative  mean  kinetic  energy  at  the  top  of  the  roughness  element. 

PkUk*/PeUe*«  1*  the  key  physical  parameter  In  evaluating  the  tendency  of 

roughness  Induced  disturbances  to  promote  boundary  layer  transition. 

Laminar  boundary  layer  profiles  are  a required  Input  to  the  evaluation  of  this  hypo- 
thesis and  they  are  also  essential  to  the  related  extension  of  the  disturbance  modeling 
to  conditions  of  Interest  fbr  flight.  For  application,  however,  the  transition  criteria 
must  be  formulated  In  terms  of  Information  available  during  the  Integral  solutions  used 
In  shape  change  codes.  This  Information  does  not  Include  boundary  layer  profiles.  The 
practical  utility  of  the  above  hypothes.1s  Is  realized  primarily  In  relating  the  effects 
of  ablation  to  the  well  understood  effects  of  variation  In  wall  cooling  ratio.  Once 
this  Is  accomplished,  the  empirical  disturbance  parameter  can  be  generalized,  In  terms 
available  during  an  Integral  solution,  to  Include  the  effects  of  ablation. 

In  the  next  subsection,  the  technique  used  to  represent  laminar  kinetic  energy 
profiles  Is  discussed.  This  Is  followed  by  an  analysis  of  these  profiles  (for  condi- 
tions typical  of  the  wind  tunnel  calorimeter  experiments)  which  substantiates  the 
hypothesis.  The  hypothesis  is  then  assumed  to  also  apply  In  the  pretence  of  blowing 
and  Is  used  In  the  determination  of  a modified  empirical  disturbance  parameter  which 
Includes  the  effect  of  blowing.  The  Influence  of  mass  transfer  on  the  boundary  layer 
profiles  Is  determined  by  the  dimensionless  blowing  rate,  B'  ■ 

Finally,  using  the  same  ''.pproach,  this  disturbance  modeling  Is  extended  to  Include 
the  effects  of  dissociation  and  ablation  typical  of  graphite  ablation  In  flight. 


A. 4.1  Boundary  Layer  Profile  Representation 


Laminar  boundary  layer  solutions  for  the  situations  of  Interest  were  generated 
using  the  Boundary  Layer  Integral  Matrix  Procedure  (BLIMP)  code  which  Is  described  In 
detail  In  Reference  A-8.  BLIMP  Is  essentially  an  Implicit  coda  for  the  solution  of 
the  two-dimensional  or  axlsymmetrlc  boundary  layer  equations  for  transport  of  mass, 
momentum,  energy,  end  species,  assuming  chemical  equilibrium.  Nonsimllar  terms  In 
the  equations  are  reduced  significantly  by  transformation  to  the  appropriate  Levy- 
Lees  coordinates. 

For  an  axlsynnetrlc  geometry  these  coordinates  are 
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In  th*  study  of  tht  solutions  dsscribid  bolow  It  Is  convenient  to  plot  boundary 
layer  profiles  In  terms  of  a normal  coordinate  which  absorbs  the  Influence  of  the 
boundary  layer  edge  conditions  but  does  not  reflect  the  Influence  of  the  surface 
boundary  cond1t1ons»  T^/T,  and  B'.  The  effects  of  varying  these  parameters  will 
appear  In  the  profiles  and  not  In  the  coordinates.  The  incompressible  form  of  the 
above  normal  variable. 
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Is  the  coordinate  which  provides  these  useful  features  and  It  Is  utilized  below  In 
the  presentation  of  boundary  layer  profiles.  It  should  be  emphasized  that  this  coor- 
dinate Is  simply  convenient  for  the  Interpretation  of  solutions  and  Is  not  Involved 
In  the  generation  of  solutions. 

For  the  purposes  of  comparing  the  relative  Influence  of  blowing  and  wall 
cooling  ratio  on  the  kinetic  energy  profile,  solutions  for  the  stagnation  region, 
where 

® " 2u“ 


are  used  In  all  that  follows.  This  Is  a good  approximation  for  the  subsonic  region 
to  which  the  present  data  and  analysis  are  restricted.  Note  that 

J^no(y.e)- 

Figure  A-13  shows  relative  kinetic  energy  profiles  for  air  In  a stagnation 
region  without  mass  transfer,  since  the  Initial  goal  Is  to  Interpret  data  taken  In 
a wind  tunnel  with  T^  « 1300*‘R,  the  enthalpy  level  for  these  solutions  was  selected 
accordingly.  Note  that  at  some  fixed  distance  from  the  wall  the  variation  of  kinetic 
energy  with  wall  temperature  Is  In  the  direction  consistent  with  trends  In  the  transi 
tion  data.  That  Is,  for  a higher  wall  temperature  ratio,  the  relative  kinetic  energy 
and,  hence,  the  magnitude  of  the  disturbance  Is  lower  so  that  transition  Is  less 
likely. 

A. 4. 2 Correlation  of  the  Empirical  and  Phsvclal  Disturbance  Parameters.  No  Blowing 

In  this  section,  the  relative  kinetic  energy  at  the  top  of  the  roughness  ele- 
ments, for  B'  ■ 0,  Is  shown  to  be  a function  of  (T./T^)(k/9)  only;  that  Is,  there  Is 
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very  little  explicit  dependence  of  p^^Vp^u^*  on  T^T^. 
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Thi  rougtinnt  •Itmnts  art  treated  at  If  th«y  protrude  a distance  k/2  Into  a 
smooth  wall  boundary  layer.  This  Is  essential 1y  the  distance  Vrom  the  mean  surface 
location.  Fbr  selected  values  of  (T^/T^^)(k/e)  and  T^/T^,  the  dimensionless  location 
of  the  top  of  the  roughnMS  elements  In  the  boundary  layer  may  be  written  as 

The  dimensionless  momentum  thickness, 


Is,  for  B'  ■ 0 and  B ■ 0.5,  a function  of  T^/T^  only.  This  function,  evaluated  from 
BLIMP  solutions  Is  shown  In  Figure  A*14.  Also  shown  are  values  of  this  function 
calculated  by  Dewey  and  Gross  (Reference  A-9)  who  report  Integral  parameters  and 
wall  derivatives  for  similar  laminar  boundary  layer  solutions.  For  laminar  boundary 
layers  with  no  blowing,  such  close  agreement  (~1  percent)  Is  expected. 

Now,  the  kinetic  energy  may  be  picked  off  the  appropriate  cooling  ratio  curve 
In  Figures  A-13.  The  results  of  this  exercise  are  shown  In  Figure  A*15.  Clearly,  the 
two  disturbance  parameters  correlate  very  well  with  each  other.  The  conclusions  Is 
that  the  transition  Reynolds  numbers  (Re^)  could  have  been  correlated  with  the  kinetic 
energy  ratio  about  as  well  as  with  the  empirical  disturbance  parameter 

(T^/T^)(k/e).  Indeed  the  correlation  shown  In  Figure  A>15  Indicates  that  the  wind 
tunnel  data  substantiate  the  hypothesis,  stated  above,  regarding  the  key  physical 
measure  of  the  disturbance. 

Other  profile  parameters,  Uj^/u^,  P|(U(j/pg«g»  and  were  examined  In 

addition  to  the  kinetic  energy  ratio  as  candidates  for  the  physical  measure  of  the 
disturbance.  None  of  these  correlates  well  with  the  empirical  measure  of  the  dis- 
turbance, (T^/T^^)(k/e),  and  therefore  all  were  rejected. 

A.4.3  Effects  of  Blowing  on  the  Roughness-Induced  Disturbance 

Because  the  kinetic  energy  profile  adequately  represents  the  effect  of  wall 
temperature  ratio  variation  on  the  rough  wall  disturbance.  It  Is  now  assumed  that  the 
modifications  to  the  kinetic  energy  profile  associated  with  blowing  would  likewise 
be  reflected  In  a modified  rough  wall  disturbance.  Therefore,  boundary  layer  calcu- 
lations Including  both  blowing  and  wall  cooling  ratio  variations  were  carried  out 
In  order  to  provide  kinetic  energy  profiles  from  which  to  assess  modified  disturbance 
magnl tudes. 


Figurts  A-16«,  b,  and  c show  tha  Influanca  of  air  to  air  blowing  on  tht  kinatic 
anargy  proflla  for  tha  wall  cooling  ratio  ranga  of  primary  Intarast.  Nott  that  blowing 
dacraasas  tha  kinatic  anargy  at  any  fixad  an  affect  which  Is  similar  to  tha  affact 
of  Incraasing  tha  wall  cooling  ratio. 

Sinca  boundary  layer  profiles  era  not  usually  generated  during  engineering 
reentry  heat  transfer  calculations,  the  relative  kinetic  anargy  Itself  Is  not  really 
suitable  for  use  directly  In  a transition  criterion.  Rather,  It  Is  more  useful  to 
reverse  the  procedure  which  lad  to  the  establishment  of  tha  relative  kinetic  energy 
ratio  as  the  physical  measure  of  the  disturbance.  The  object  then  Is  to  establish 
a correlation  between  the  physical  disturbance  parameter,  P|(U|(*/pgU.^,  and  some 
modified  empirical  disturbance  parameter,  say 

1 k 

Since  replaces  the  wall  cooling  ratio  In  the  empirical  disturbance  parameter,  It 
must  reduce  to  the  cooling  ratio  for  B*  ■ 0. 

The  function  <F(T^/T^,B‘)  may  be  defined  as  that  value  of  (TM/T0)gt.Q  which 
leads  to  the  tame  relative  kinetic  energy  at  the  same  position  In  the  boundary  layer. 
For  several  discrete  values  of  ngi  lines  of  constant  B*  were  plotted  on  a relative 
kinetic  energy  versus  wall  cooling  ratio  plane.  One  of  these  plots  Is  shown  In 
Figure  A>17  where  some  examples  of  ip  determinations  are  Illustrated.  Values  of  i|i 
from  several  such  plots  are  shown  correlated  In  Figure  A>18.  The  straight  lines 
shown  have  the  equation 


Figure  A-19  substantiates  this  result.  The  correlation  obtained  for  no  blow- 
ing In  Figure  A-1S  remains  valid  for  the  solutions  which  Include  blowing  when  the 
empirical  disturbance  parameter  Is  modified  as  Indicated  above.  The  relative  kinetic 
energy  for  each  point  shown  on  the  graph  Is  determined  from  tha  plots  In  Figures  A-16a, 
b,  and  c,  using  the  no  blowing  dimensionless  momentum  thickness  evaluated  at  the 
actual  wall  cooling  ratio  to  calculate  ng(k/2). 

A.4.4  Implications  for  Graphite  Ablation  In  Flltfit 

The  primary  differences  between  the  flight  and  the  low  temperature  air  to  air 
blowing  situations  Just  discussed  are  the  differences  In  chemical  composition  arising 
from  dissociation  and  graphite  ablation.  With  molKular  weight  varying  through  the 
boundary  layer,  the  density  distribution  Is  no  longer  the  Inverse  of  the  temperature 
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Figure  A-19.  Cerreletlm  of  tke  relative  kinetic  energy  vltli  the 
■otfiflod  e^lrlcal  dlstnvhance  panaeter. 
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distribution.  Sinct  donslty  apposrs  dirtctly  In  tht  kinotic  tntrgy.  It  would  statn 
moro  appropriate  to  uto  the  density  ratio.  P,/p^.  In  the  empirical  disturbance  param- 
eter rather  than  the  wall  cooling  ratio. 

1Wo  addi  ^nal  BLIMP  calculations  were  performed  In  order  to  assess  this 
proposed  p^/p^  i>  dificatlon.  A graphite  ablation  solution  was  generated  for  a 
typical  reentry  situation  at  50.000  feet  assuming  chemical  equilibrium  and  a steady 
state  surface  energy  balance.  Using  the  wall  temperature  found  for  that  case,  a non- 
ablating solution  was  generated  so  that  the  Influence  of  dissociation  could  be 
examined  Independent  of  ablation.  These  solutions  were  used  to  test  the  correlation 
between  relative  kinetic  energy  and  the  disturbance  parameter  (k/t|)8)  us1,ng  alternate 
choices  for  4*.  Figure  A-20  shows  these  results  compared  to  the  correlation  line 
from  Figure  A- IS. 

For  the  nonablating  solution,  the  correlation  using  the  wall  cooling  ratio  Is 
unsatisfactory  except  near  the  outer  edge  of  the  boundary  layer.  Using  the  density 
ratio  the  points  near  the  surface  are  brought  nicely  into  line  with  the  correlation 
obtained  for  low  enthalpy  conditions. 

It  also  appears  that  the  density  ratio  properly  correlates  the  dimensionless 
momentum  thickness.  For  this  case,  p^/p^  ■ 0.439,  while  ng(6)  ■ 0.431.  This  Is  In 
agreement  with  the  results  presented  In  Figure  A-14  when  the  abscissa  Is  Interpreted 
as  the  density  ratio  (T^/T^  « p^/p^  for  the  solutions  on  which  this  figure  Is  based). 

When  the  density  ratio  Is  used  In  conjunction  with  the  blowing  term  developed 
earlier,  the  ablating  solution  Is  also  well  correlated.  Although  this  technlqli'  does 
not  yield  a good  correlation  near  the  edge  of  the  boundary  layer,  this  region  Is  not 
too  Important  because,  for  flight  conditions  of  Interest,  roughness  heights  are  less 
than  subsonic  boundary  layer  thicknesses.  Thus,  the  best  choice  for  the  disturbance 
Is  k/(4i6),  where 


The  density  ratio,  as  used  hare.  Includes  the  effects  of  ablation.  This  quantity 
Is  not  readily  available  from  engineering  codes  but  Its  calculation  Is  within  the 
scope  of  this  type  of  code.  For  the  case  of  graphite  ablation,  use  of  the  air 
density  at  the  correct  wall  temperature  would  result  In  k/ij;e  being  about  10  percent 
low.  This  In  turn  results  In  a 2 kft  lower  transition  altitude  prediction. 
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Figure  A-20.  OHmlatlon  of  tfw  rclotivo  kinetic  eeergy  with  the  wodlflcd  dis 
porawetcr  for  f11^  cowHtlows. 


In  SMiMryt  transition  location  and  onsat  crltarla,  obtained  as  correlations 
of  solid  calorimeter  wind  tunnel  data,  have  been  rationally  extended  to  Include  the 
flight  case  of  carbon  ablation  at  high  enthalpy.  Accurate  prediction  of  transition 
during  flight  should  now  be  possible  provided  sufficiently  accurate  evaluation  of 
the  surface  roughness  height  Is  performed.  Since  fabrication  variations  apparently 
produce  significant  differences  In  surface  roughness  within  a particular  type  of 
graphite  (ATJ-S.  for  example),  the  determination  of  roughness  height  should  be  car- 
ried out  using  a similarly  prepared  sample.  Furthermore,  the  appropirate  roughness 
Is  the  one  produced  after  soam  amount  of  laminar  ablation. 
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